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ABSTRACT 

We identify and investigate the nature of the 20 brightest 250 /im sources detected by 
the Balloon-borne Large Aperture Submillimetre Telescope (BLAST) within the central 
ISOarcmin^ of the GOODS-South field. Aided by the available deep VLA L4GHz radio 
imaging, reaching 51.4 — 40 /iJy (4-cr), we have identified radio counterparts for 17/20 of 
the 250 fim sources. The resulting enhanced positional accuracy of 1 arcsec has then al- 
lowed us to exploit the deep optical (HST), near-infrared (VLT) and mid-infrared (Spitzer) 
imaging of GOODS-South to establish secure galaxy counterparts for the 17 radio-identified 
sources, and plausible galaxy candidates for the 3 radio-unidentified sources. Confusion is a 
serious issue for this deep BLAST 250 /im survey, due to the large size of the beam. Never- 
theless, we argue that our chosen counterparts are significant, and often dominant contributors 
to the measured BLAST flux densities. For all of these 20 galaxies we have been able to de- 
termine spectroscopic (8) or photometric (12) redshifts. The result is the first near-complete 
redshift distribution for a deep 250 /im- selected galaxy sample. This reveals that 250 /im sur- 
veys reaching detection limits of 40mJy have a median redshift z 1, and contain not 
only low-redshift spirals/LIRGs, but also the extreme z 2± 2 dust-enshrouded starburst galax- 
ies previously discovered at sub-millimetre wavelengths. Inspection of the LAB OCA 870 /im 
imaging of GOODS-South yields detections of c:^ 1/3 of the proposed BLAST sources (all at 
z > 1.5), and reveals 250/870 /im flux-density ratios consistent with a standard 40 K modified 
black-body fit with a dust emissivity index (3 = 1.5. Based on their IRAC colours, we find 
that virtually all of the BLAST galaxy identifications appear better described as analogues of 
the M82 starburst galaxy, or Sc star-forming discs rather than highly obscured ULIRGs. This 
is perhaps as expected at low redshift, where the 250 /im BLAST selection function is biased 
towards spectral energy distributions which peak longward of A^es^ = 100 /im. However, it 
also appears largely true dX z ^2. 
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1 INTRODUCTION 

The observed far-infrared background peaks around a wavelength 
A 200 - 250 /im (Puget et al. 1996; Fixsen et al. 1998). 
Over the past 10 years, deep extragalactic surveys have closed in 
on this wavelength regime from both shorter and longer wave- 
lengths, and there has been dramatic progress in uncovering pop- 
ulations of sources which undoubtedly contribute to this back- 
ground. Specifically, at wavelengths ^ 4 times longer than the 
peak (A ^ 850 — 1200 /xm), surveys with the Sub-millimetre 
Common-User Bolometer Array (SCUBA), Mambo, AzTEC and 
now LABOCA (e.g. Coppin et al., 2006; Bertoldi et al. 2007; Weiss 
et al. 2009; Austermann et al. 2010) have been used to success- 
fully assemble substantial samples of submm- selected extragalac- 
tic sources, and the spectral energy distributions (SEDs) of these 
sources are certainly rising steeply towards shorter wavelengths. 
Conversely, at wavelengths ~ 4 times shorter than the peak, the 
Spitzer Space Telescope has proved effective at extragalactic source 
selection up to wavelengths A ^ 70 /xm (Magnelli et al., 2009). But 
until now, the effective production of deep extragalactic samples in 
the central wavelength range A :^ 100 — 500 /xm has been prohib- 
ited by the atmosphere and/or limitations on instrument sensitivity 
and telescope aperture. 

This situation should shortly be transformed by the SPIRE in- 
strument on the 3.5-m diameter Herschel Space Observatory (Grif- 
fin et al. 2007). However, a powerful first insight into the nature of 
the 250 /im-selected galaxy population is already being provided 
by the the Balloon-borne Large Aperture Submillimetre Telescope 
(BLAST). 

BLAST is a 1.8 m diameter stratospheric balloon telescope 
that operates at an altitude of approximately 35 km, above most 
of the atmospheric water vapour which essentially prohibits effec- 
tive far-infrared observations from the ground. BLAST is thus half 
the size of Herschel, and it is equipped with a camera which is a 
prototype of the SPIRE camera, enabling simultaneous broad-band 
imaging with central wavelengths of 250, 350 and 500 /xm. In 2006, 
BLAST undertook an 11 -day flight from Antarctica, during which 
it executed an observing programme which included the first deep 
far-infrared imaging survey ever undertaken within the Extended 
Chandra Deep Field South (ECDFS). 

This unique survey (described in more detail in Section 2) 
has already been the subject of several investigations. First results 
on the far-infrared number counts and the resolution of the back- 
ground were presented by Devlin et al. (2009), with more detailed 
investigations of these two key topics being presented by Patan- 
chon et al. (2009) and Marsden et al. (2009). Pascale et al. (2009) 
used the survey to constrain cosmic star-formation history, while 
Viero et al. (2009) studied correlations in the background. Dye et 
al. (2009) and Ivison et al. (2010) used the supporting radio (VLA) 
and mid-infrared Spitzer imaging of the field to attempt to iden- 
tify countparts to the brighter BLAST sources, and to explore the 
far-infrared/radio correlation in distant galaxies. 

The study by Dye et al. (2009) showed that the BLAST beam, 
at least at 250 /xm, is still small enough to allow the successful 
identification of a reasonable fraction of BLAST sources with radio 
and mid-infrared counterparts. However, the supporting data over 
the full 9 deg^ BLAST survey area are of variable quality, and Dye 
et al. attempted to secure identifications for sources selected at all 
3 BLAST wavelengths. As a result, they succeeded in identifying 
counterparts for only c::^ 55% of the sources in the BLAST cata- 
logues detected at > 5 cr in at least one waveband. This inevitably 
limits the conclusions that can be drawn about the redshift distri- 



bution of the BLAST source population, although Dye et al. con- 
cluded that 75% of the sources in the shallow+deep catalogue 
lie at z < 1. 

The aim of the study presented here is to establish a near- 
complete redshift distribution for a subset of the BLAST 250 /im 
sources selected down to a fainter flux-density limit of ^so — 
40 mJy, and subsequently to establish the basic physical properties 
(stellar mass, size, morphology, star-formation history, AGN con- 
tent) of the galaxies which host the observed 250 /xm emission. To 
try to achieve this, we have deliberately confined our attention to 
the small sub-area of the BLAST map which contains the very best 
optical {HST ACS), near-infrared (VLT ISAAC) and mid-infrared 
{Spitzer) imaging data, and the highest density of optical spectro- 
scopic redshifts. This area, the 150 arcmin^ GOODS-South field 
(Dickinson et al. 2004), is also where the radio (VLA) imaging is 
most sensitive, and it has been mapped in its entirety at 870 /xm by 
LABOCA (Weiss et al. 2009), and at 1.1 mm by AzTEC (Scott et 
al. 2010). GOODS-South is also the chosen location for some of 
the first deep high-resolution near-infrared imaging currently be- 
ing undertaken with WFC3 on the refurbished HST, and is the fu- 
ture target of planned deep sub-mm and far-infrared imaging with 
SCUBA2 (Holland et al. 2006) and SPIRE-hPACS on Herschel 

Thus, by confining our attention to GOODS-South we can ex- 
plore just what can be achieved given the best possible support- 
ing data. Conversely, we can also establish what level of support- 
ing data is actually required for an effective, complete study of 
the future, larger, 250 /xm-selected galaxy samples which will be 
produced by Herschel. In addition, the wealth of existing informa- 
tion on other galaxy populations in this field (see Section 2) facil- 
itates the comparison of BLAST galaxies with reference samples 
of field galaxies selected at similar redshifts and/or stellar masses, 
and to explore the frequency of mergers. Finally, the recent acqui- 
sition of deep sub-millimetre and millimetre wavelength imaging 
in the same field allows a first exploration of how the 250 /im and 
870/1100 /xm populations are related. 

To assemble a useful sample of 250 /im sources within this 
field, we have pushed the BLAST source selection significance 
threshold to ~ 3.5cr. However, the resulting sample is still suffi- 
ciently small (20 sources) that we can afford to explore the prop- 
erties of individual sources in some detail. As a result this work 
should be regarded as a pilot study of the mix of source popula- 
tions which can be uncovered in a deep 250 /im survey, and an 
exploration of how best to overcome the problems that are encoun- 
tered in trying to identify and study the faint far-infrared sources in 
maps which are inevitably highly confused (note that the Herschel 
beam at 500 /im is essentially the same size as the BLAST beam 
at 250 /im). The broader statistical robustness of the results pre- 
sented here is clearly limited both by small number statistics, and 
by cosmic variance. 

The layout of the paper is as follows. In Section 2 we de- 
scribe the selection and robustness of the BLAST 250 /im sub- 
sample within the central region of GOODS-South, and then sum- 
marize the wealth of available multi-frequency imaging and spec- 
troscopy in the field. In Section 3 we describe the process of 
obtaining radio and optical counterparts for the BLAST sources, 
and highlight the care required to overcome the problems associ- 
ated with source blending in the confused BLAST maps. In Sec- 
tion 4 we provide or derive redshifts for all of the sources, in- 
cluding secondary/alternative counterparts. Then, in Section 5, we 
present multi-colour images of the BLAST galaxies and provide 
notes on each of the individual sources. Finally, in Section 6 we 
review our results and discuss the robustness of our conclusions. 
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In particular we present and discuss the implications of the strik- 
ing associations found between the BLAST 250 jim. sources and 
the LABOCA 870 iim. sources recently uncovered by Weiss et al. 
(2009). We also discuss the robustness of our derived redshift distri- 
bution for 250 /im-selected sources, and briefly compare the mid- 
infrared colours of BLAST galaxies with those of other known 
galaxies at both high and low redshifts. Our conclusions are sum- 
marized in Section 7. 



2 DATA 

2.1 The BLAST 250 yun source sub-sample 

The 250 jiui source sample selected for study here represents a 
small subset of the far-infrared (250 /xm, 350 /xm, 500 fim) sam- 
ple uncovered in the ECDFS by the BLAST Antarctic flight of 
December 2006 (Devlin et al. 2009). An area of 8.7 deg^ centred 
on RA03^ 32"", Dec - 28° 12' (J2000) was mapped by BLAST 
to mean sensitivities of (7250 = 36mJy, (7350 = 30mJy, and 
(7500 = 20 mJy per beam. Within this 'BLAST GOODS South 
Wide' (BGS-Wide) map, a deeper image, reaching (7250 = 11 mJy, 
cTsso = 9 mJy, and (7500 = 6 mJy per beam was made of the cen- 
tral 0.8 deg^ centred on the southern field of the Great Observa- 
tories Origins Deep Survey (GOODS; Dickinson et al. 2003) at 
RA03^ 32™ 30", Dec - 28° 48' (J2000). It is from the core of 
this 'BLAST GOODS South Deep' (BGS-Deep) map that we have 
extracted a 'complete' sub-sample of the BLAST 250 fim. sources 
which lie within the sub-region with the very best supporting multi- 
frequency data. As shown in Fig. 1, we have chosen to define this 
area as the ^150 sq. arcmin region of sky covered by the GOODS- 
South HST ACS 4-band optical imaging (Giavalisco et al. 2004). 

Because our primary aim here is to study the optical/infrared 
galaxy counterparts of the BLAST sources, we have confined our 
attention to 250 /im sources, because the larger BLAST beam sizes 
at 350 and 500 /im make it even more difficult to reliably iden- 
tify all but the brightest longer-wavelength sources (the BLAST06 
beams are best fit by Gaussians with FWHM of 36, 42, and 
60 arcsec at 250 /xm, 350 /xm and 500 /xm respectively; Pascale et 
al. 2008). Our study thus differs from, and complements that of 
Dye et al. (2009), who sought identifications for all 351 BLAST 
sources detected at > 5cr in any one of the three BLAST wavebands 
across the whole of BLAST GOODS South survey. We have also 
extracted a significantly deeper sub- sample of sources for study, 
selecting all 250 fim sources which lie within the region shown in 
Fig. 1 and have a signal: noise ratio greater than 3.25. Pushing the 
BLAST data in this way allows us to construct a sample of useful 
size from within the relatively small area covered by the highest 
quality supporting datasets, and also allows a first investigation of 
the nature of the fainter 250 /im source population which will be 
studied in detail by Herschel. 

Imposing this signalinoise cut produced a sample of 24 
sources. We chose this significance threshold because, below this 
level, there are only 9 sources with signal: noise > 2, and only three 
of these were reproduced at > 2 cr in an alternative map reduction 
(and all three at < 2.8cr). In other words, this threshold appears 
to mark a clear (not unexpected) point below which source relia- 
bility degenerates markedly, at least in this particular region of the 
BLAST map. 

Of these 24 sources, we then rejected three (176, 228, 777) as 
inappropriate for this study because their nominal positions moved 
by greater than 20 arcsec between different map reductions. This 



does not necessarily mean these sources are not real peaks in the 
map (they range in S/N from 5cr to 3.5cr), but simply that their 
positions appear to be too sensitive to the precise choice of map- 
production/source-extraction parameters to allow reliable determi- 
nation of the most likely radio or optical counterpart (possibly be- 
cause they are extreme examples of confusion from a number of 
fainter sources). There certainly appears to be something different 
about these sources because all other sources in our sample moved 
by less than 5 arcsec between alternative map reductions. Finally, 
we also rejected the second least significant source (1158), because 
it was not found at > 2 cr in the alternative map reduction. 

This leaves a 'clean' sample of 20 sources for further study. 
The 250 /xm positions of the BLAST sources are indicated in Fig. 1, 
and tabulated in Table 1, along with their 250 //m flux densities 
and signalinoise ratios. Given that the 4 rejected sources displayed 
a range of 250//m signalinoise, the remaining sample can be re- 
garded as representative down to a flux-density detection limit of 
6*250 — 35 mJy in the zero-mean map (5*250 — 45 mJy relative 
to the local background - see below). The mean signalinoise ratio 
for these 20 sources is 5.6, but this is biased by the brightest 2 or 3 
sources. The median signalinoise is 4.2, and this number is adopted 
as representative where appropriate later in this paper. 

Before searching for identifications, it is worth pausing to con- 
sider the robustness of the far-infrared sources listed in Table 1. 

The 250 /im flux densities given in Table 1 are the 'raw' val- 
ues as extracted from the BLAST map via convolution with the 
relevant noise- weighted point- spread function (Truch et al 2009), 
effectively a weighted fit of the beam to the raw unsmoothed data 
to derive the best (maximum likelihood) estimate of the flux density 
of a point source in the map (e.g. Serjeant et al. 2003). Due to the 
steepness of the 250 /xm source counts, and the large BLAST beam, 
the flux densities of the individual sources (which would be seen in 
a higher resolution map) will undoubtedly have been boosted by the 
effects of confusion (i.e. flux contributions from fainter sources) 
and by Eddington bias (a combination of noise and steep source 
counts), especially for the lower signalinoise detections (Scott et 
al. 2002; Coppin et al. 2006; Pantachon et al. 2009). Reading fluxes 
from the zero-mean map can, however, partially offset the flux con- 
tributions from fainter confused sources, which in effect constitute 
the local background for any brighter source in the map. However, 
detailed simulations are clearly required to establish robust statisti- 
cal corrections for sources at a given flux density (Patanchon et al. 
2009). 

The other issue to consider is the noise. The flux-density noise 
levels quoted in Table 1 include only the instrumental noise which, 
within this deep section of the map, is typically cr25o — llmJy 
per beam. However, analysis of the pixel flux-density distribution 
within the relevant 150 arcmin^ region under study here, shows that 
the spatial noise in the map is (7250 — 15mJy per beam (as mea- 
sured from a Gaussian fit to the negative tail of the flux-density dis- 
tribution). This implies that the confusion noise is comparable to 
the instrumental noise, at cr25o ^10 mJy per beam. The situation 
is therefore similar to that encountered with the deepest 850 /im 
maps made of the Hubble Deep Field with SCUBA on the JCMT 
(Hughes et al. 1998; Peacock et al. 2000). 

It might therefore seem appropriate to increase the noise level 
to cr25o — 15 mJy per beam when assessing the significance of our 
sources in the actual map. However, if we choose to do this, for 
consistency we must also adopt the local background in the map 
which examination shows to be typically —11 mJy per beam. This 
is not surprising, as the mode of the distribution of flux-density 
values in the zero-mean map lies at —12 mJy per beam. Thus, per- 
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Figure 1. The location of the 20 BLAST 250 /im sources in GOODS-South superimposed on the H ST ACS z-bmd imaging mosaic of the field (LH plot) and 
on a greyscale of the BLAST 250 /im map itself (RH plot). Sources are numbered by BLAST source ID (from the full BLAST map), as given in Table 1, with 
higher numbers indicating lower 250 /im flux densities. The circles have a radius of 15 arcsec, which is the radius adopted for the search for radio and optical 
counterparts to the BLAST sources. The actual BLAST beam has a FWHM of 36 arcsec at 250 /im 



forming aperture photometry on the map typically increases flux- 
densities by 11 mJy per beam, and boosts the assumed noise level 
to (7250 — 15mJy per beam. Adoption of these two corrections 
reduces only slightly the formal significance level of the sources 
listed in Table 1, and still leaves our faintest source above the 3a 
threshold. 



Given these complications, in this paper we make little di- 
rect use of the 250 /im flux densities, other than to discuss the 
use of 250/870 /im flux-density ratio as a consistency check on de- 
rived source redshift (see Section 6.1). We also do not attempt to 
exploit the 350 or 500 /im measurements of our 250 /im- selected 
sources, as the larger beams, and relatively low sensitivity (at least 
to 250 /im-selected sources) of BLAST at these longer wavelengths 
prohibits the extraction of usefully accurate flux densities. We thus 
do not quote 350 or 500 /im flux densities for the sources in Table 
1 (although values for the brighter sources can be found in Dye 
et al. 2009). We do, however, check, on a source-by- source ba- 
sis, for detections or non-detections at these longer wavelengths, 
as described in Section 5. Finally, we note that as there are ^ 500 
BLAST 250 /im beams in the 150 arcmin^ area under study here, 
we would expect 0.5 non-existent sources to clear our adopted 
3.3cr threshold purely by chance. Relative to the 250 /im map mode 
of — 12mJy, the MUSIC region under study here contains one 3a 
negative peak. We thus conclude that, while flux boosting may be 
substantial for individual sources, our adopted signaLnoise thresh- 
old is reasonable. 



2.2 Supporting multi-frequency data 

2.2.1 Radio: VLA 1.4-GHz imaging 

As demonstrated by the follow-up of sub-millimetre surveys with 
SCUBA, very deep VLA imaging is necessary (and frequently suf- 
ficient) for the successful identification of the galaxy counterparts 
of sources detected with the large beams of current far-infrared/sub- 
millimetre facilities (e.g. Ivison et al. 2002, Ivison et al. 2007). This 
works for three reasons. First, star-forming galaxies produce copi- 
ous quantities of synchrotron emission. Second, even in the deepest 
available radio maps, 1.4 GHz sources still have sufficiently low 
surface-density that positional coincidences within a 'reasonable' 
search radius are generally statistically rare. Third, if a robust radio 
counterpart is found, the 1 arcsec positional accuracy provided 
by the VLA at 1.4 GHz essentially always yields an unambiguous 
optical/infrared galaxy counterpart for further study. 

Very deep {aiA — 7.5 /iJy), high-resolution 1.4 GHz imag- 
ing of GOODS-South is now available at the centre of the ECDFS 
as described by Kellermann et al. (2008) and Miller et al. (2008). 
As discussed by Dye et al. (2009), the published source catalogue 
adopts a very conservative detection threshold of la. We therefore 
re-analysed the image of Miller et al. (2008) using the techniques 
described by Ibar et al. (2009) to create a radio catalogue down to a 
4 a limit of 30 /iJy (at which depth the cumulative source density 
on the sky is ^ 0.8 arcmin"^). This catalogue was then searched 
for potential radio counterparts to the BLAST sources using the 
method described below in Section 3. 
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BLAST SOURCE NAME BLAST RA250 DEC250 S/N S'250 

(lAU) ID (J2000) (J2000) (250) /mJy 



BLAST J033235- 


-275530 


(250 


/im) 


4 


53.146220 


-27.925272 


16.2 


177 


BLAST J033229- 


-274414 


(250 


/im) 


6 


53.123900 


-27.737326 


14.4 


157 


BLAST J033218- 


-275216 


(250 


/im) 


59 


53.075139 


-27.871343 


6.7 


74 


BLAST J033205- 


-274645 


(250 


ptm) 


66 


53.022647 


-27.779397 


6.6 


72 


BLAST J033235- 


-274932 


(250 


ptm) 


104 


53.147540 


-27.825642 


6.0 


65 


BLAST J033217- 


-275054 


(250 


/im) 


109 


53.072586 


-27.848349 


5.9 


64 


BLAST J033221- 


-275630 


(250 


/im) 


158 


53.088715 


-27.941911 


5.2 


58 


BLAST J033212- 


-274642 


(250 


/im) 


193 


53.051434 


-27.778541 


4.9 


54 


BLAST J033217- 


-274944 


(250 


/im) 


257 


53.073511 


-27.829046 


4.6 


51 


BLAST J033242- 


-275514 


(250 


/im) 


318 


53.176619 


-27.920633 


4.4 


48 


BLAST J033238- 


-274954 


(250 


/im) 


503 


53.159514 


-27.831733 


4.0 


44 


BLAST J033213- 


-274302 


(250 


/y,Tn^ 


552 


53.054825 


-27.717303 


3.9 


43 


BLAST J033247- 


-275418 


(250 


/im) 


593 


53.197956 


-27.905028 


3.9 


43 


BLAST J033243- 


-275146 


(250 


/im) 


637 


53.181931 


-27.862863 


3.8 


42 


BLAST J033232- 


-275304 


(250 


/im) 


654 


53.136852 


-27.884531 


3.8 


42 


BLAST J033243- 


-274650 


(250 


/im) 


732 


53.179377 


-27.780768 


3.7 


40 


BLAST J033213- 


-274246 


(250 


/im) 


830 


53.056853 


-27.712811 


3.6 


40 


BLAST J033248- 


-274443 


(250 


/im) 


861 


53.200852 


-27.745501 


3.6 


39 


BLAST J033232- 


-274558 


(250 


/im) 


983 


53.134915 


-27.766112 


3.4 


38 


BLAST J033246- 


-275321 


(250 


/im) 


1293 


53.193438 


-27.889410 


3.3 


36 



Table 1. The BLAST 250 /im-selected sample in GOODS-South. The positions are derived from optimal beam fitting to the BLAST map. 6*250 is the flux 
density derived from the zero-mean BLAST map. 



In cases where the radio emission appears heavily resolved 
(e.g. for the brightest source, BLAST 4, where at first sight there 
appear to be 6 individual VLA candidate counterparts) we mea- 
sured the total radio flux density using TVS TAT within AIPS. We 
checked these values against the flux densities measured with the 
16-arcsec beam of the Australia Telescope Compact Array (ACTA) 
(Norris et al. 2006), and found them to be in good agreement for 
the 4 sources which proved sufficiently bright to have also been 
detected at 1.4 GHz by ACTA (Afonso et al. 2006). 

2.2.2 Optical: HST ACS imaging 

Deep optical imaging over the 150 arcmin^ area shown in Fig. 1 has 
been obtained with the Advanced Camera for Surveys (ACS) on 
board HST in 4 different filters: F435Ty (^435), ^6061^ (Veoe), 
F77bW (^775) and F850LP (zgbo)- These data were taken as part 
of GOODS (Giavalisco et al. 2004) and reach 5 a Hmiting (AB) 
magnitudes (within a 1 arcsec diameter aperture) of 26.9, 26.9, 
26.1, and 25.8, respectively. By definition, all 20 of the BLAST 
sources lie within the area covered by the GOODS-South HST 
ACS optical imaging. One of the BLAST sources (BLAST 732) 
lies within the even deeper ACS imaging provided by the Hub- 
ble Ultra Deep Field (HUDF: Beckwith et al. 2006), over an area 
of 11 arcmin^ centred on RA03^ 32"" 39s, Dec - 27° 47' 29.1'' 
(J2000). 

2. 2. 3 Near- infrared data: VLT and HST imaging 

Deep near-infrared (J, if, K) imaging of almost all of the HST ACS 
field illustrated in Fig. 1 has now been completed with the ISAAC 
camera on the VLT (Retzlaff et al. 2010). This imaging covers 
^ 143 arcmin^, of which ^ 136 arcmin^ overlaps with the HST 
ACS imaging. All but 2 of the BLAST sources under study here lie 
within the area covered by the ISAAC images. 

This near-infrared imaging is of excellent quality, with the 
FWHM of the PSF having a median value of ^ 0.5 arcsec (it varies 



from ~ 0.37 to ^ 0.7 arcsec across the field; Bouwens et al. 2008). 
Partly as a result of this good image quality, the point source sensi- 
tivity of this imaging is very deep, reaching 5 a detection levels of 
J 25, and H, Ks ^ 24.2 (AB magnitudes). 

As described by Bouwens et al. (2008), various HST NIC- 
MOS Jiio and ifieo imaging programs have been undertaken in 
GOODS-South. However, the coverage offered by this imaging 
is too patchy to be of much use in the present study, and frus- 
tratingly (but not unexpectedly) no BLAST source lies within the 
5.8 arcmin^ ultra-deep NICMOS sub-image of the HUDF under- 
taken by Thompson et al. (2005). However, BLAST 732 does lie 
within the new Wide Field Camera 3 (WFC3) lios, -^125, i^ieo 
imaging of the HUDF, which reaches a 5-cr limiting magnitude 
(AB) of 29 in all three near-infrared bands (McLure et al. 2010; 
Oesch et al. 2010; Bouwens et al. 2010). 



2.2.4 Mid-infrared: Spitzer imaging 

Again as part of GOODS, ultra-deep Spitzer imaging with the In- 
frared Array Camera (IRAC: Fazio et al. 2004) has been obtained 
over the whole of the ACS GOODS-South field illustrated in Fig. 1 
(proposal ID 194, Dickinson et al., in preparation), in all 4 IRAC 
channels (3.6, 4.5, 5.6, and 8.0 /xm). The IRAC 5 a detection limits 
are S3.6 ^ 25.9, S4.5 ^ 25.5, S5.6 ^ 23.3, Sg.o - 22.9 (AB 
magnitudes). 

The 24 /xm Spitzer MIPS data originally obtained as part of 
GOODS has been augmented and incorporated within the Spitzer 
Far-Infrared Deep Extragalactic Legacy (FIDELfl survey (Mag- 
nelli et al. 2009), and reach 5-cr detection limits of S24 — 30 /xJy. 



PI M. Dickinson, see http : / /vrvrvr . noao . edu/noao/f idel/] 
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2.2.5 Sub-mm: LABOCA 870 fim. survey 

The full 30 X 30 arcmin ECDFS has recently been mapped at a 
wavelength of 870 /im by the LABOCA ECDFS Sub-miUimetre 
Survey (LESS) (Weiss et al. 2009). The 12-m diameter of the 
APEX telescope (Gusten et al. 2006) delivers a 19.2-arcsec 
FWHM beam. The LESS image has a uniform depth of as 70 — 
1.2 m Jy beam" ^ (as measured in the map, and hence includ- 
ing confusion noise) and the positions and flux densities of the 
LABOCA sources have been determined in a similar manner to 
that described above for the BLAST source extraction (Weiss et al. 
2009). 

The combined analysis of the BLAST and LABOCA LESS 
source lists and maps in GOODS-South is the subject of a major 
study which will be presented elsewhere (Chapin et al. 2010). Here 
we have confined our attention to checking which of the 20 BLAST 
250 /im sources in Fig. 1 are also present in the LESS 870 /im cat- 
alogue presented by Weiss et al. (2009), and briefly exploring the 
implications of the detections and non-detections (see Section 6.1). 

2.2.6 Spectroscopic redshifts 

In recent years the GOODS-South field has been the target of a 
series of spectroscopic campaigns (Dickinson et al. 2004, Stanway 
et al. 2004, Strolger et al. 2004, van der Wei et al. 2004 Szokoly et 
al. 2004, Le Fevre et al. 2004, 2005, MignoH et al. 2005, Doherty et 
al. 2005, Roche et al. 2006, Ravikumar et al. 2007, Vanzella et al. 
2009, Tayler et al. 2009, Eales et al. 2009). Of particular potential 
importance for the present study has been the completion of the 
ESO/GOODS VLT/FORS2 programme in the GOODS-South field 
(Vanzella et al. 2005, 2006, 2008). 

The combined impact of these programs is that over 1000 
galaxies in the GOODS-South field shown in Fig. 1 now possess 
reliable spectroscopic redshifts. As a result, and unsurprisingly, it 
transpires that all except one of the BLAST galaxy counterparts 
which we identify below with a low-redshift (i.e. z < 1) galaxy 
already possesses a known spectroscopic redshift. 

However, what is perhaps surprising, especially given the 
work of Vanzella et al., is that only one of the potential higher- 
redshift BLAST galaxy counterparts currently has a spectroscopic 
redshift (BLAST 1293, z = 1.382; Vanzella et al. 2008). As ex- 
plored further below, this is primarily a consequence of the fact 
that the high-redshift BLAST galaxy counterparts are too faint/red 
for optical redshifts to be determined, even in the VLT/FORS2 
campaign. The exploitation of the aforementioned multi- waveband 
photometry for the production of photometric redshift estimates 
thus remains important for this and future work on the study of 
sources selected at far-infrared/sub-millimetre wavelengths. 

2.2.7 Photometric redshifts 

Several photometric redshift catalogues have now been published 
for the ECDFS in general and for the central GOODS-South field 
in particular. 

The COMBO 17 project (Wolf et al. 2004, 2008) covers the 
whole 30 X 30 arcmin area of the ECDFS, and exploits 17-band (5 
wide and 12 narrow) optical imaging obtained with the Wide Field 
Imager (WFI) at the ESO 2.2-m telescope. As explained by Wolf 
et al. (2004, 2008) these 'very low resolution spectra' allow solid 
{cFz — 0.02 — 0.05) photometric redshifts to be obtained down to 
a magnitude Hmit of R 23. While the COMBO 17 catalogue 
contains estimated redshifts reaching z :^ 2, in practice the depth 



of the optical imaging, and lack of near-infrared information limits 
its usefulness to ^ < 1 where the optical spectroscopy in GOODS- 
South is very complete. Consequently, a COMBO 17 redshift was 
only adopted for one BLAST source (BLAST 257-1, z = 0.689). 

Over the 30 x 30 arcmin ECDFS, photometric redshift cat- 
alogues based on the available broad-band optical+near-infrared 
data have been produced by the MUSYC (Multiwavelength Sur- 
vey by Yale-Chile) collaboration (e.g. Taylor et al. 2009), and (for 
24 /xm sources) by Rowan-Robinson et al. (2008). However, of 
more importance for the current study are the photometric redshift 
catalogues assembled for the deep central 150 arcmin^ area by 
the GOODS-MUSIC project (Grazian et al. 2006), by the alter- 
native analyses of Caputi et al. (2006) and Dunlop, Cirasuolo & 
McLure (2007), and most recently by Brammer et al. (2008), ap- 
plying the EAZY algorithm to the FIREWORKS data (Wuyts et 
al. 2008). 

We utilised results from these studies but, as described in Sec- 
tion 4, we also decided to derive new redshift estimates for each of 
the BLAST source candidate identifications. This was partly neces- 
sary simply because not all of the sources feature in existing photo- 
z catalogues. However, it also allowed us to check the impact of 
deriving photometry through different apertures, and to explore the 
effect on Zphot and Szphot of allowing extreme values of extinction 
Av (see Dunlop, Cirasuolo & McLure 2007). 



3 SOURCE IDENTIFICATION 
3.1 Radio identifications 

There is no doubt that, even with the smallest BLAST beam, 
(FWHM=36 arcsec at 250 /im) identifying secure unambiguous 
optical/IR galaxy counterparts to the far-infrared sources is a chal- 
lenge. It is also clearly the case that some (possibly large) subset of 
the BLAST sources may be the product of the blended far-infrared 
flux density from two or more galaxies at the same, or differing 
redshifts. 

For some of the brighter BLAST sources identified by Dye 
et al. (2009) and Ivison et al. (2010), sufficiently bright and rare 
24 //m counterparts exist to allow the unambiguous identification 
of the BLAST sources with Spitzer mid-infrared sources. How- 
ever, at the fainter flux densities which are the subject of the 
present study, the sheer number-density of MIPS sources in the 
deep Spitzer imaging makes the isolation of statistically significant 
unambiguous mid-infrared counterparts virtually impossible. 

We have therefore focussed here on the search for associations 
with VLA 1.4 GHz sources which, even at the aforementioned deep 
flux-density limit of 5'i.4 :^ 30 /iJy, remain sufficiently rare on the 
sky(A^(5'i.4 > 30/iJy) 0.8 arcmin"^) for random associations 
to be statistically unlikely. 



3.1.1 Selection of candidate radio counterparts 

We searched for radio- source counterparts within a radius of 
15 arcsec of each BLAST 250 pm. source. This choice of search 
radius is not arbitrary, and can be justified in a number of ways. 

First, we note that the anticipated positional uncertainty in 
a source detected with an (undeboosted) signahnoise ratio, S/N, 
with a Gaussian beam of FWHM 0, in the presence of background 
power-law cumulative number counts of the form N(> f) oc 
is given by (Ivison et al. 2007): 
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Table 2. Galaxy identifications derived from associations with 1.4 GHz VLA sources (17 objects) or analysis of the optical-infrared photometry of possible 
candidates (3 sources) within a search radius of 15 arcsec. 6 sources have more than one statistically acceptable galaxy counterpart (given in parenthesis). Note 
that = refers to the sum of the individual values of P for only the 17 primary (i.e. lower P value) radio identifications. 
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Into this expression we insert the FWHM of the BLAST 250 //m 
beam ^ = 36 arcsec, our modal value of S/N = 4.2, and an assumed 
faint count power-law index 7 = 1 (slightly sub-Euclidean num- 
ber counts appear appropriate at fainter flux densities, as a double 
power-law fit transits from 7 = 2.5 to 7 = 0.8 at ^so — 30mJy; 



Patanchon et al. 2009). The result is cTp 



6 arcsec. 15 arcsec is 



then the 2.5cr search radius which is expected to contain ^ 95% 
of all genuine radio counterparts (see Ivison et al. 2007). Note here 
that the number count correction term is applied to effectively de- 
boost the source flux density for Eddington bias, and has the effect 
of reducing a raw S/N = 4.2 to an effective deboosted S/N = 3.4. 

Second, while expanding the search radius further might result 
in one or two additional radio counterparts, the increased number 
of ways in which a counterpart can be selected inevitably weakens 
the statistical security of associations which could still have been 
found within a smaller search radius (see discussion of probabil- 
ity of mis-identification below). From a series of simulations we 
have established that, given the depth of the radio data available in 
this field, the number of statistically secure radio identifications is 
optimised by adoption of a search radius between 12 and 15 arc- 
sec. Fortunately this number is (just) consistent with the 2 to 2.5 a 
positional uncertainty in the BLAST 250 jim. sources under con- 
sideration here. This provides some retrospective justification for 
our decision to confine our study of the faint sources to the BLAST 



250 /im catalogue. It would not be possible to successfully under- 
take this type of analysis for ^ 4cr sources if the beam was much 
bigger than 36 arcsec (as is the case, for example, at 500 //m). 

Third, after determining the most likely galaxy identifications, 
one can check that the distribution of angular offsets is consistent 
with that expected given the anticipated positional uncertainty in 
the BLAST sources. We perform this check below, after careful 
consideration of the confusing effects of source blending. 

Finally, we experimented with scaling the search radius with 
BLAST source S/N, but did not find that this significantly affected 
the results. In addition, such refinement of the search radius may 
place too much confidence in the precise accuracy of the positions 
of the brighter BLAST sources. This is not to say there is any evi- 
dence for serious systematic errors in the BLAST positions. Indeed 
the positional offset of only 1.5 arcsec to the galaxy counterpart of 
the brightest BLAST source in the field provides confidence that 
any additional pointing errors in the BLAST 250 fim. positions are 
minimal. Rather, this is dangerous because, as discussed below (see 
also Ivison et al. 2010), many of the brightest sources may consist 
of blends. In this situation, blind application of equation (1) can 
result in statistically robust primary and secondary radio identifica- 
tions being missed because of an inadequate search radius (e.g. for 
BLAST 6, rigorous application of equation (1) would yield a l.Scr 
search radius of 3.8 arcsec, and a failure to find either of the two 
contributing galaxies - see Section 3.1.3). 
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Figure 2. Left: Cumulative distribution of angular offsets between BLAST 250 and radio positions for the 17 radio identifications, compared with the expected 
distribution for 18 sources given the typical deboosted S/N of the BLAST detections and the size of the 250 /im Gaussian beam (i.e. a = 6 arcsec). Right: The 
equivalent plot after moving the radio position to the mean of the radio position of the two alternative IDs for the 4 BLAST sources 6, 59, 503 and 830, which 
all appear to be potential 2-source blends within the BLAST beam. Moving these 4 sources restores the expected behaviour, showing that the distribution of 
offsets is as expected. 



3.1.2 Calculation of probability of mis-identification 

Having chosen a candidate radio identification within our 15 arcsec 
search radius, we then calculate the probability that such a coin- 
cidence could have occurred by chance. Following Downes et al. 
(1986; see also Dunlop et al. 1989), we calculated the raw Pois- 
son probability that a radio source of the observed 1.4 GHz flux 
density would be discovered by chance at the measured distance 
from the nominal BLAST source position. We then correct this a 
priori probability for the number of ways such a statistical coin- 
cidence could have been discovered given the available search pa- 
rameter space defined by the maximum search radius (15 arcsec), 
the radio- source number density at the limiting search flux-density 
(N{SiA > 30 /iJy) ~ 0.8 arcmin"^), and the form of the radio 
source counts over the flux-density range of interest (here we adopt 
a power-law index of 1.4). For the deep data under study here, this 
correction is often substantial, typically increasing P by a factor of 
a few. This technique has been applied previously to estimate the 
robustness of radio identifications for SCUBA sources (e.g. Ivison 
et al. 2002, 2007), and should yield similar results to the method 
adopted by Dye et al. (2009) based on Monte Carlo simulations. 

The value of P thus calculated is the probability that the ob- 
served association is the result of chance. We stress that even a very 
low value of P does not prove that the radio source is the BLAST 
source. Nevertheless, a low value of P clearly does imply that the 
radio source is likely related to the BLAST source in some way. 
This could be true for several reasons. First, the radio source could 
be the single true counterpart of the BLAST source. Alternatively 
it could be one of group of 2 or more galaxies which contribute to 
the BLAST flux-density peak (of particular relevance here, given 
the large beam-size). Third, the statistical result could be a con- 
sequence of some secondary association (e.g. clustering with the 
true BLAST sources, or the result of gravitational lensing). Despite 



these worries over statistical interpretation, we re-emphasize that 
our search for radio counterparts is not motivated purely by their 
statistical rarity and good positional accuracy, but also by the phys- 
ical evidence that dust-enshrouded star-forming galaxies also pro- 
duce copious quantities of synchrotron emission (resulting in the 
well-known far-infrared:radio luminosity correlation - e.g. Ivison 
et al. 2010). A sensible hypothesis, therefore, is that an apparently 
statistically secure radio counterpart to a BLAST source is at least 
a significant contributor to the observed 250 //m emission. 



3.1.3 Results 

The results of this process are tabulated in Table 2. Radio counter- 
parts were found for 17 of the 20 BLAST sources, with multiple 
counterparts found for 5 of these 17. Half of the BLAST sources 
have radio counterparts with P < 0.05, and all 17 radio counter- 
parts have P < 0.1. While it would be nice to insist on a signifi- 
cance threshold of P < 0.05, we here adopt P < 0.1 as the best 
that can be realistically achieved given the large search radius re- 
quired by the BLAST beam. Reassuringly, a sum of the P values 
for the 17 primary identifications yields only SP c± 0.8, suggest- 
ing that ^ 1 source has been misidentified. At the same time our 
adopted search radius would lead us to expect to have missed ^ 1 
true radio identification at larger radii, so our failure to find a radio 
counterpart for 3 BLAST sources is not altogether surprising. The 
technique used to isolate the possible optical counterparts for these 
3 remaining BLAST sources (listed in Table 2) is explained below 
in subsection 3.2 (although we note here the possibility that 1 or 2 
of these 250 /im of these sources could be erroneous - see Section 
6). 

As a consistency check, in the left-hand panel of Fig. 2, we 
plot the cumulative distribution of angular offsets between BLAST 
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250 /im and radio/optical positions for the 17 primary identifica- 
tions listed in Table 2, and compare this with the 'expected' distri- 
bution for the afore-mentioned circular Gaussian with a = 6 arcsec 
(assuming 18 sources in total, to allow for the fact our search ra- 
dius should only include 95% of identifications). While application 
of the Kolmogorov-Smirnov (KS) test reveals these distributions 
to only differ at slightly less than the 2a level, this comparison has 
clearly not worked as well as for SCUBA sources (Pope et al. 2006; 
Ivison et al. 2007). Closer inspection of Table 2 reveals that, de- 
spite low values of P, the positional offsets for some of the brighter 
BLAST+radio associations are surprisingly large. This point is per- 
haps best illustrated by comparison of the identifications for the 
brightest two BLAST sources in our sub-sample. For BLAST 4, 
the radio identification is completely secure (P = 0.0002) and, as 
expected, lies only 1.5 arcsec from the BLAST 250 /im position 
(predicted positional uncertainty apos = 1.35 arcsec for this 16- 
cr BLAST source). By contrast, the most likely identification for 
BLAST 6, while still statistically compelling (P = 0.004), Hes 11 
arcsec from the BLAST 250 fim position, a 7 cr deviation in posi- 
tional offset for a 14 cr BLAST source. 

The key to resolving this apparent discrepancy lies in the fact 
that BLAST 6 has a second (also statistically compelling) radio 
identification lying almost equi-distant from the 250 /im position, 
in the diametrically opposite direction. It turns out that BLAST 
6 consists of an interacting pair of galaxies, whose positions (as 
marked by their VLA centroids) are separated by 20 arcsec, and 
which both lie ai z = 0.076 (see Section 4). Clearly the radio iden- 
tification procedure adopted here has yielded the correct redshift, 
but the position of the 250 /im source has either been distorted by 
contributions from both galaxies, or arises from dust which lies in 
the region between the two optical galaxies. Either way, the ap- 
parently excessive distance to either alternative radio identification 
should not be naively interpreted as a failure of the radio identifi- 
cation technique, or as casting doubt on the BLAST positional ac- 
curacy. Thus, in the case of multiple identifications, each of which 
appears statistically secure, a fairer assessment of the identifica- 
tion procedure is provided by comparing the BLAST position with 
the mean position of the two alternative radio identifications, es- 
pecially when such identifications straddle the nominal 250 //m 
position. This is in fact the case for 4 of the 5 BLAST sources 
listed in Table 2 which have two alternative radio identifications - 
BLAST 6, 59, 503 and 830. If, for these 4 sources, we replace the 
primary identification radio source position with the mean radio 
position of both alternative radio counterparts, the cumulative dis- 
tribution of positional offsets changes to that presented in the right- 
hand panel of Fig. 2. This shows that allowance for the possibil- 
ity of 250 /im source blending in just these 4 apparently multiply- 
identified sources, is sufficient to bring the distribution of angular 
offsets completely into line with expectation. 

This analysis provides increased confidence in the statistical 
robustness of our identification procedure. However, it does serve 
to highlight the potential problems of confusion arising from the 
large size of the BLAST beam. In particular it shows that, even 
when our identification procedure can lead to an unambiguous red- 
shift, it is unclear what fraction of the 250 /xm emission can be reli- 
ably attributed to a given radio identification. To explore this further 
we attempted to distinguish whether, even for a clear-cut high S/N 
case such as BLAST 6, we could distinguish whether or not the 
BLAST source was a single source, or a blend of both galaxies. 
The result, summarized and discussed in Fig. 3, is that we cannot, 
even when we are confident that the positions of the potentially 
blended sources are well known, and separated by more than half 



the FWHM of the BLAST beam. Thus it is clear that we need to be 
very cautious in attempting to combine 250 /xm and 1.4 GHz data 
to determine, for example, the far-infrared:radio flux ratio for indi- 
vidual BLAST sources, despite the fact that for only two sources in 
the current sample (BLAST 59 and BLAST 503) does the choice 
between two alternative radio identifications significantly influence 
the inferred redshift (see Section 4). There are clear lessons here 
that care will need to be taken in deriving the far-infrared spectral 
energy distribution of sources uncovered by Herschel, especially at 
500 /im where the beam-size is comparable to that which applies 
here at 250 /im. 

Finally, in Fig. 4 we plot the radio flux-density distribution 
of the primary identifications listed in Table 2. This figure demon- 
strates that most of the radio identifications have Si a < 100 /iJy 
and that several have flux densities close to the 30 /iJy limit. This is 
another reason that the existence of 3 radio-unidentified sources in 
our 20-source sample should not be regarded as surprising. Indeed, 
Fig. 4 makes the generic point that near-complete radio identifi- 
cation of a 250/im-selected galaxy sample requires 1.4 GHz radio 
data reaching at least the same sensitivity in /i Jy as is achieved at 
250//m in mJy. 

3.2 Optical/infrared galaxy counterparts 

For the radio-identified sources, determining the correct 
optical/near-infrared galaxy counterpart is then relatively straight- 
forward, due to the high accuracy of the VLA positions. Most, 
but not all, of these galaxies are fisted in the GOODS-MUSIC 
catalogue, and we give GOODS -MUSIC ID numbers for each 
object in Table 3. For those sources which do not feature in 
GOODS MUSIC we were able to successfully identify the galaxy 
counterpart in the i^^s-band and IRAC imaging (one source, 593, 
is completely undetected in the ACS optical imaging, even at 
2:850 ), and then performed our own photometry using the available 
multi-frequency imaging. 

This leaves the 3 radio-unidentified BLAST sources, 193, 654, 
and 732. We must bear in mind the real possibility that, given the 
S/N ratio of the source catalogue under study, at least one of these 
250 /j,m sources may not be real. However, as discussed in more de- 
tail in the source notes (see Section 5), 193 appears to be confirmed 
at 350 /im, 654 appears to be confirmed by a detection at 500 /im, 
and 732 has apparently been detected at 870 /im by LABOCA. We 
also checked that modest extension of the radio search radius (to 
20 arcsec) would still not have yielded a possible radio counter- 
part for these sources. Thus it seems more likely that these sources 
may have radio detections lying just below the current radio image 
threshold, possibly because they lie within the high-redshift tail of 
the source population. 

We therefore attempted to establish 'best-bet' optical/infrared 
counterparts for these sources by deriving a photometric redshift 
for every possible galaxy counterpart within the same 15 arcsec 
search radius used for the radio identifications, and then confin- 
ing our attention to candidates with z > 1.5. The reason for re- 
stricting the potential redshift range to z > 1.5 is that all three 
of these BLAST sources, if real, have 250 //m/1.4GHz flux ratios 
^25o/*S'i.4 > 1000, a value which is not exceeded by any known 
galaxy spectral energy distribution until z > 1.5. As described in 
more detail in section 5, this yields unique candidates for 193 and 
654, and 3 candidates for 732. The positions and GOODS-MUSIC 
ID numbers (where available) for these objects are also listed in 
Tables 2 and 3. 

The identifications for these sources must inevitably be re- 
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Figure 3. Results of an attempt to deblend the 250 jim image of BLAST 6. 
At the top we show a 200 x 200 arcsec B+V+z HST ACS colour image of 
the field centred on BLAST 6. The upper 2 grey-scale panels then both show 
the same region as imaged by BLAST at 250 /xm (with 10 arcsec pixels). 
The single cross in the upper-left image shows the position of a proposed 
single source which when convolved with the BLAST beam and fitted to 
the data yields a best-fit flux density of 5*250 = 157 mJy. When subtracted 
from the data this leaves the residual image shown in the lower-left panel. 
The two crosses in the upper-right image show the proposed positions of 
two distinct sources (i.e. the positions of the radio sources 6-1 and 6-2) 
which, when convolved with the BLAST beam and fitted to the data, yield 
best-fit flux densities of S'250 = 91mJy (lower source, 6-1 in Table 2) 
and S'250 = 75mJy (upper source, 6-2 in Table 2). When subtracted from 
the data this combined-source model leaves the residual image shown in 
the lower-right panel. Unfortunately, due to the BLAST beam, even in this 
rather well-defined test case it is not possible to distinguish whether the 
single or double source model is a better description of the BLAST data, 
as the single-source fit yields = 520, and the double-source fit yields 
= 530 (for 400 pixels). Thus, while the redshift of this source is not in 
doubt, and it is clear that both galaxies are related to the 250 /im emission, 
it is not possible to decide whether this emission arises primarly between 
the two optical galaxies, or is a blend of emission from both nuclei. 



Sl.4GHz//^Jy 

Figure 4. The distribution of 1 .4 GHz radio flux density, in 20 /iJy-wide 
bins, for the radio-identified sources in the 20-source BLAST 250 /im 
GOODS-South sample (excluding the two very bright sources at z < 0.1). 
All radio flux densities are as given in Table 2. The dotted line marks the 
4cr 30 /iJy limit of our radio candidate search in the VLA data. Given this 
distribution of flux densities, and expecially the fact that several identifi- 
cations have flux densities just above the 30 /iJy limit, the existence of 3 
radio-unidentified sources in our 20-source sample is not surprising. This 
figure demonstrates that near-complete radio identification of a 250 12m- 
selected galaxy sample requires 1.4 GHz radio data reaching at least the 
same sensitivity in /i Jy as is achieved at 250 fim in mJy. 

garded as more speculative than most of the radio identifications 
(indeed we later reject 654 on the basis of no detection at 870/im 
or 24/im). However, the selected objects provide a plausible galaxy 
counterpart for each of these sources which is at least consistent 
with the other available photometric and redshift infomation. 

The optical, near-infrared and Spitzer photometric data for all 
27 galaxies listed in Tables 2 and 3 are presented in Appendix A. 



4 REDSHIFTS 

4.1 Spectroscopic redshifts 

A search of the latest redshift catalogues resulting from the nu- 
merous spectroscopic surveys within GOODS-South (see Section 
2.2.6) produced redshifts for 13 of the radio-identified galaxy coun- 
terparts listed in Table 2. These redshifts are presented in Table 3, 
where it can be seen that 8 are for primary counterparts, and 5 are 
for secondary counterparts. The highest spectroscopic redshift is 
z 1.382. 

4.2 Photometric redshifts 

For 637-2 we adopted the COMBO 17 redshift estimate of 
z = 0.689, but no other galaxy counterpart lacking a spectro- 
scopic redshift has a robust COMBO 17 redshift. Reference to the 
GOODS -MUSIC catalogue revealed that this is because the re- 
maining candidates all appear to lie at 2; > 1. 
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Rather than simply adopting the redshift estimates for the 
higher redshift galaxies from the existing published catalogues, we 
chose to re-analyse the photometry of each proposed identification 
to derive our own values of Zphot- We chose to do this for 5 reasons. 
First, not all of the candidates are listed in the existing catalogues, 
because they are too faint, and for these we had to extract our own 
photometry from the imaging data and perform the first estimate of 
their redshifts. Second, given we are potentially dealing with very 
dusty galaxies, we wished to explore the effect on Zphot of allow- 
ing extinction (Av)to range to large values (as previously explored 
in a different context by Dunlop et al. 2007). Third, we wished to 
to marginalise over the full parameter space to derive meaningful 
confidence intervals on the estimated redshifts (errors on individual 
redshift determinations are frequently not provided in the published 
catalogues). Fourth, we simply wanted to check the robustness of 
the published redshift estimates, and explore the effect of extracting 
the photometry through, for example, apertures of different sizes. 
Fifth, we wish to derive physical parameters such as galaxy age and 
stellar mass, to facilitate the further study of the physical properties 
of the BLAST galaxies. 

The SED fitting procedure we applied to derive the photomet- 
ric redshifts and physical properties of the BLAST galaxy counter- 
parts is based largely on the public package Hyperz (Bolzonella 
et al. 2000). The observed photometry was fitted with synthetic 
galaxy templates generated with the stellar population models of 
Chariot & Bmzual (e.g. Bmzual 2007). We used a variety of star- 
formation histories: instantaneous burst and exponentially declin- 
ing star-formation with e-folding times 0.1 < T(Gyr) < 10, as- 
suming solar metallicity and a Salpeter initial mass function (IMF). 
For dust reddening we adopted the prescription from Calzetti et al. 
(2000) within the range < Ay < 4. We also included absorption 
due to HI clouds along the line of sight in the intergalactic medium, 
according to Madau (1995). 

Our best estimate of Zphot for every source (including those 
with spectroscopic redshifts) is tabulated in Table 3. For every 
galaxy which lacks a spectroscopic redshift we show in Fig. 5 
the best-fitting spectral energy distribution and a plot of ver- 
sus redshift marginalised over all other fitted parameters (age, star- 
formation history, dust extinction, and normalization). The la con- 
fidence interval in redshift is given below each plot. 

As a test of the robustness of our photometric redshifts we 
compare our values for Zphot with the available spectroscopic red- 
shifts in Fig. 6. This only provides a direct test out to z 1.5, 
but confirms good agreement and reveals no catastropic outliers. In 
this same figure we also compare our results for the higher redshift 
sources with the results in the GOODS -MUSIC catalogue, where 
such values are available. While the scatter in the comparison is 
inevitably larger, the results still provide confidence that the values 
of Zphot for the higher redshift sources are also reasonably robust. 

In Fig. 7 we provide multi-wavelength postage stamp images 
of the primary BLAST identifications, ranked by redshift. Similar 
visual information is provided for the secondary/alternative coun- 
terparts in Fig. 8. 

The final redshift distribution for the GOODS-South 250 //m 
'brightest-counterpart' sample is presented in Fig. 10, and dis- 
cussed in Section 6.2 

5 NOTES ON INDIVIDUAL SOURCES 
BLAST 4: Zspec = 0.038 

The brightest 250 /xm source in the field is unambiguously 
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Table 3. Spectroscopic and photometric redshifts for the BLAST galaxy 
identifications (including alternative identifications listed in parenthesis). 
Only for 2 of the 6 sources with alternative identifications does the choice 
significantly affect the inferred redshift (59 and 503). Also given is the 
sub-millimetre flux density (aSsto) or 3-a limit taken from the LABOCA 
LESS survey of GOODS-South by Weiss et al. (2009). Of the 11 proposed 
BLAST sources with z > 1.5, LABOCA has detected 6. In contrast, no 
870 /im detections were achieved for any of the BLAST sources identified 
with galaxies at z < 1.5. As discussed in Section 6, this is largely as ex- 
pected, and provides an independent sanity check on the basic form of the 
redshift distribution presented here. 

identified with the brightest optical/near-infrared galaxy in 
GOODS-South. This is a very low-redshift, edge-on disc galaxy 
in which the effect of patchy dust obscuration can clearly be seen 
in the HST ACS colour composite image shown in the first panel 
of Fig. 7. Comparison with the VLA imaging apparently reveals 
5 statistically secure alternative radio counterparts (VLA sources 
499, 524, 572, 701, 986, all of which have P < 0.05), but in fact 
all of these lie near to the centroid of the optical galaxy, and it is 
clear that the VLA imaging has fragmented a more extended radio 
source. We have adopted VLA source 499 as the formal radio ID 
(position listed in Table 2); this is both the closest candidate and 
has the smallest value of P (and is the same radio ID adopted by 
Dye et al. 2009). A comparison of the sum of the flux densities 
of the 5 VLA counterparts {^Sia = 202 mJy) with the flux 
density derived from the lower resolution (16 arsec FWHM) ACTA 
imaging of the field (ACTA Sia = 310 =b 48 fiJy; Norris et al. 
2006) confirms that much of the radio emission from the galaxy 
may have been resolved out by the VLA imaging. Our own sum of 
all the flux-density in the VLA map yields 5'i.4 = 680 =b 34 /iJy, 
and we adopt this as the best available estimate of 1.4 GHz flux 
density. This lowers the inferred S'250 '■ Sia flux density ratio to 
^ 260 from the erroneously high value of 3000 implied by the 
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numbers quoted in Dye et al. (2009). Given the angular size of 
the galaxy, and the lack of deep low-resolution radio imaging with 
a beam size comparable to BLAST, this ratio is still probably an 
over-estimate of the true value. The fact that the BLAST position 
lies only 1.5arcsec from the centroid of the optical/radio galaxy 
provides confidence in the BLAST 250 /im positional accuracy 
because, at S/N =16, the rms positional error for a BLAST 250 //m 
source is expected to be 1.35 arcsec (see Section 3.1). Thus, at least 
in this region of the 250 /xm map, any systematic pointing error 
appears impressively small. The source is also clearly detected at 
350 /im and 500 /im, although the 500 /im position lies 15 arcsec 
distant from the true galaxy centroid. This is not hugely unexpected 
given the rms positional error of 7.5 arcsec anticipated from a 5.5cr 
detection with the 60-arcsec FWHM 500 //m BLAST beam, but 
serves as a reminder of why it is virtually impossible to securely 
identify BLAST 500 /im sources which lack 250 //m counterparts. 

BLAST 6: z^pec = 0.076 

The second brightest 250 )Ltm source in the field is clearly associ- 
ated with an interacting galaxy pair. Both the 250 yL/m and 350 //m 
positions appear to lie almost equi-distant between two compa- 
rably bright, apparently interacting galaxies at the same redshift 
(see Fig. 1 and Fig. 3). The result is two statistically significant 
alternative radio IDs (VLA sources 72 and 100), and we have 
adopted the southern source (VLA 72) as the chosen ID because 
it has a marginally smaller value of P (the 500 /xm position also 
favours the southern galaxy). It is this galaxy which is shown in 
Fig. 7, with the alternative galaxy identification shown in Fig. 8. 
Because both objects lie at the same redshift, the precise choice 
obviously does not actually affect our final redshift distribution at 
all. However, in truth the inferred positional offset of 10.4 arcsec 
between the radio and 250 /xm positions is unexpectedly large for 
a 14 cr BLAST source and a 40cr VLA source. Thus as discussed 
in detail in Section 3.1.3, it is essentially certain that the 250 /im 
emission either receives comparable contributions from both 
galaxies, or in fact does actually arise from a region between the 
two optical galaxies. Unfortunately, as shown in Fig. 3, the BLAST 
data do not allow us to distinguish between these two scenarios. 
As with BLAST 4, the integrated 1.4 GHz flux density of this 
source (quoted in Table 2 for both galaxies) is much larger than the 
nominal value derived from the VLA radio catalogue by Dye et al. 
(2009). 

BLAST 59: Zest = 2.29 (or Zspec = 1.097) 
Stepping down by a factor of two in 250 /im brightness, the third 
brightest BLAST source has a flux density of S'250 — 75mJy, 
more typical of what is expected in a field of this size given 
the source counts (Devlin et al. 2009). However, for this object, 
selection of the correct identification is extremely difficult. The 
BLAST centroid lies roughly equidistant between 3 alternative 
radio counterparts. The statistically preferred radio ID for this 
source lies to the east of the 250 /xm position, but gains additional 
support from the 500 /im position. Statistically, the association 
seems reasonably secure. The radio source has an optical/infrared 
counterpart in the GOODS MUSIC catalogue (shown in Fig. 7), but 
this has no spectroscopic redshift. The GOODS MUSIC estimated 
redshift is z = 2.13, while our own analysis of the photometry 
yields z — 2.29 (2.05 — 2.35) as shown in the first panel of Fig. 5. 
However, given the position of the 250 /im source it is extremely 
likely that at least some of the far-infrared flux is contributed by 
the other radio-identified objects, which are galaxies lying within 
a known large-scale structure at z :^ 1.09. We thus also retain (as 



59-2) the second most likely radio identification, VLA 240, which 
is a galaxy ai z = 1.097 which displays AGN emission, and lies to 
the southwest of the 250 //m centroid (Fig. 8). This is the one case 
in the sample where both alternative radio counterparts also appear 
to be associated with individual LAB OCA sources, strengthening 
the case that the 250 /xm emission results from the blending of at 
least these two sources. 

BLAST 66: Zest = 1 94 

This 250 /im source has clear detections at 350 /im and 500 //m, 
both reassuringly close to the 250 /im position. The radio identifi- 
cation is robust and unambiguous. The source lies just outside the 
ISAAC near-infrared imaging, but the ACS and Spitzer IRAC pho- 
tometry is adequate to yield a unique and well-constrained redshift 
solution at 2: c:^ 2, as shown in Fig. 5. The high-redshift nature 
of this source receives support from its detection by LAB OCA 
at 870 /xm. The multi-colour postage stamps in Fig. 7 reveal a 
faint complex source at optical wavelengths, which is bright and 
unconfused at 3.6 /im. This is the third of the three sources in 
this deep GOODS-South sample which were also identified by 
Dye et al. (2009), and not surprisingly they deduced the same 
VLA source as the correct ID. Confusingly however, they list a 
C0MB017 redshift of z = 1.16 which is clearly inconsistent 
with the photometric redshift constraints derived here. However, 
upon inspection of the latest COMBO 17 catalogue it can be seen 
that the estimated redshift of this source (C0MB017 35066) is 
highly uncertain (unsurprisingly so, given that the object is very 
faint, with R 25.8) and that the peak of the redshift probability 
distribution is listed as 2; = 1.84, consistent with our value. The 
redshift listed by Dye et al. (2009) can thus be safely rejected. 

BLAST 104: ^spec = 0.547 

This 250 //m source has detections at 350 /im and 500 /im which 
reinforce the BLAST position. The VLA identification is centred 
on a red galaxy with a spectroscopic redshift of 2; = 0.547. The 
closest companion seen in Fig. 7 lies at Zest = 1.44, but the 
object ENE of the BLAST position has Zest = 0.54 and lies 
within 28kpc of the BLAST source at this redshift (although it 
has no radio detection, so presumably is contributing less to the 
far-infrared flux density). The next nearest radio counterpart lies 
15.7 arcsec north-east of the BLAST position, and thus just outside 
our adopted search radius. This object has a spectroscopic redshift 
z = 2.578. 

BLAST 109: Zspec = 0.124 

This 250 /im source is not found in the BLAST catalogues at 
longer wavelengths, suggesting it lies at only moderate redshift. 
There are two formally- significant radio counterparts, but one 
lies close to our search boundary at d = 15 arcsec. We adopt the 
closer radio counterpart, with the lower value of P as the most 
likely identification. This ties the BLAST source to a low-redshift, 
edge-on disc galaxy. 

BLAST 158: Zest 1 85 

Secure 350 /xm and 500 /xm counterparts to the 250 //m source 
confirm its position as just off the southern edge of the ISAAC 
K-hand mosaic. The radio ID is faint, but unique and reasonably 
secure. As with source 66, the available ACS optical and Spitzer 
IRAC photometry are sufficient to provide a reasonably solid 
estimated redshift for the associated galaxy, at Zest = 1-85, and 
the high-redshift nature of this source receives support from its 
detection by LAB OCA at 870 /im. In Fig. 7 the galaxy looks like 
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BLAST 59-1: Zest = 2.29 (2.05 - 2.35) 



BLAST 66: Zest = 1.94 (1.85 - 2.10) 





BLAST 158: z^st = 1.85 (1.70 - 1.95) 



BLAST 193: Zest = 1.81 (1.55 - 2.00) 



Figure 5. Photometric redshift determination for the 12 galaxy identifications in the 20-source BLAST 250 /im GOODS-South sub-sample which lack 
spectroscopic redshifts. For each source the upper plot shows the best galaxy spectral energy distribution fit to the HST optical, ISAAC near-infrared and 
Spitzer 3.6 /im and 4.5 /im photometry. The lower plot shows how varies with redshift, marginalised over galaxy age, star-formation history, and dust 
reddening (allowing the extinction to range up to Ay = 4). The BLAST ID number of each source, and its estimated redshift (with la error range) are given 
under each two-panel plot. 



an extremely complex faint system at optical wavelengths, but it is 
bright at IRAC wavelengths. 

BLAST 193: Zest = 1 81 

This 250 /im source has a solid detection at 350 /xm, but has no 
candidate VLA counterpart within our adopted search radius of 
15 arcsec (in fact none within 20 arcsec). We therefore fitted galaxy 
models to the optical-infrared photometry of all objects in the 
ISAAC K-bmid image within 15 arcsec of the 250 /xm position. 
This search produced five potential counterparts. Only one of these 
lies at z > 1.5, and we found this object to also be extremely 
red. This transpires to be one of only 6 ultra-obscured iC- selected 
galaxies in GOODS-South studied in detail by Dunlop et al. (2007; 
object 1865), and lies only 7 arcsec from the far-infrared position. 

BLAST 257: Zest = 0.689 



This 250 /im source is not found in the BLAST catalogues at 
longer wavelengths, suggesting it lies at only moderate redshift. 
This is the second source in the sample which has two alternative 
radio identifications associated with galaxies at the same redshift. 
The preferred VLA counterpart (585) has a C0MB017 redshift 
Zphot = 0.689, while the alternative radio identification (VLA 
source 852, 19 arcsec distant (so outside our formal search radius) 
lies in a galaxy with Zspec = 0.664. It seems likely that both 
galaxies contribute to the 250 /xm flux density, but the precise 
choice of ID does not affect the final redshift distribution. We adopt 
the galaxy associated with VLA source 585 as the statistically 
most likely association. 

BLAST 318: ^est = 2.09 

This 250 /im source has clear detections at 350 /im and 500 /im, 
both reassuringly close to the 250 /im position. The radio iden- 
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BLAST 552: z^st = 1.68 (1.45 - 1.80) 



BLAST 593: Zest > 2.5 (2cr) 



Figure 5. continued. 



tification is robust and unambiguous, and associated with an 
interacting galaxy pair (Fig. 7) which the photometry constrains to 
He at 2; ^ 2 (Fig. 5). The high-redshift nature of this source gains 
support from its detection at 870 /xm by LABOCA. VLT F0RS2 
spectroscopy has been attempted for this object, but yielded no 
redshift. 



BLAST 503: z^st = 1.96 (or Zspec = 0.241) 
This source is not found in the BLAST catalogues at longer 
wavelengths, which would suggests it lies at only moderate 
redshift. However, while there is a possible radio counterpart with 
a spectroscopic redshift z = 0.241 (VLA 315), the statistically 
preferred option is VLA 132. The optical to infrared images 
shown in Fig. 7 appear to reveal a complex multiple interacting 
system, which again is well constrained by the photometry to 
lie at 2; ^ 2 (Fig.7). The 250 /xm flux density may receive 
contributions from both sources, but in this case our choice of 
ID obviously does (strongly) influence the adopted redshift. 
There is no LABOCA detection to help support a high redshift, but 



neither does the non-detection strongly exclude it (see Section 6.1). 
BLAST 552: Zest = 1.68 

This source may be detected at 350 /im and 500 /im, but the 
positional agreement is poor. We have adopted VLA 932 as the 
most likely radio identification, but have concerns here about the 
extent to which the 250 /xm emission in this region is confused by 
the close proximity of BLAST 830. This is one of the more dubious 
250 /im sources in the sample. For this reason, and because this 
source does not appear to have been detected at 870 /im (see 
Section 6) it has, in the end, been excluded from the final proposed 
redshift distribution for the sample presented in Fig. 10. 

BLAST 593: ^est > 2.5 

This source has no catalogued 350 /im counterpart, but is a clear 
detection at 500 /im, and is nearly coincident with the second 
brightest LABOCA source in the field. In addition it has a solid and 
unambiguous radio counterpart in VLA 405. The host galaxy of 
this radio emission is completely invisible in the optical imaging. 
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BLAST 654: Zest = 2.62 (2.40 - 2.70) BLAST 732-1: Zest = 2.97 (2.85 - 3.05) 
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BLAST 861: Zest = 1.95 (1.62 - 2.18) BLAST 1293: Zest = 1.37 (1.35 - 1.55) 

Figure 5. continued. 

perhaps just visible in K, but clearly seen at 3.6 /xm (Fig. 7). All BLAST 654: ^^est = 2.62 

evidence points towards a high-redshift dusty galaxy (perhaps the This source has no catalogued 350 /im source, but is detected at 

most distant in our sample), but the lack of detections over a wide 500 /im only ^ 4 arcsec distant from the 250 /im position. This 

range in wavelength means that the photometric redshift is poorly 500 /im detection, combined with the lack of any radio detection 

constrained (see Fig. 5), and thus we adopt a lower limit for Zest- (no counterparts within 30 arcsec), suggest the source lies at high 

redshift. We therefore fitted galaxy models to the optical-infrared 
photometry of all objects in the ISAAC i^-band image within 

BLAST 637: Zspec = 0.279 15 arcsec of the 250 /am position. This search produced seven 

This 250 /xm source is not found in the BLAST catalogues at potential counterparts. Only one of these lies at z > 1.5, and 

longer wavelengths, suggesting it lies at only moderate redshift. it is the second closest to the 250 /im position (4.1 arcsec). It is 

There are two alternative formally- significant radio counterparts, also favoured by the 500 /im position. We have therefore adopted 

and we select the closer and marginally more significant VLA this galaxy (GOODS MUSIC 2977, at Zest — 2.62) as the best 

110 in favour of VLA 74. Both radio sources are associated with candidate identification in the available data. The GOODS -MUSIC 

galaxies with a spectroscopic redshift z = 0.279, and Fig. 7 redshift for this galaxy is z = 2.7, in excellent agreement with our 

indicates that this is another interacting galaxy pair. However, own results. The COMBO 17 redshift is z — 0.115, revealing the 

unlike the situation for BLAST 6, the far-infrared position does not severe limitations of the COMBO 17 catalogue at faint magnitudes 

lie between the two objects, and favours association with the more {R > 24). As described in Section 6, however, this source does 

luminous galaxy. Clearly the choice of galaxy counterpart does not not have a catalogued LAB OCA 870/im counterpart (as might 

influence the adopted redshift. be expect at high redshift), and it is the only primary BLAST ID 
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Figure 6. A test of the robustness of the new optical-infrared photometric 
redshifts derived here for the galaxy identifications of the 250 /im GOODS- 
South sample. Our new photometric redshifts are compared with the precise 
spectroscopic values for 8 out of the 9 sources at 2; < 1.5 (solid symbols), 
and against the pubhshed GOODS-MUSIC photometric redshifts for the 8 
other sources in the sample for which these are available (open symbols). 
The former comparison is a true test of the accuracy of our redshift esti- 
mates (at least at low redshift), while the latter comparison demonstrates 
the extent to which independent photometric redshifts agree for the fainter 
higher-redshift sources in the sample. The agreement is very good, which is 
perhaps not surprising given the generally well-defined and unique minima 
in for the individual fits shown in Fig. 5. 



with no flux at 24/xm. For these reasons, while we present the 
information we have gathered on this potential galaxy counterpart, 
we exclude it from the final redshift distribution for the sample. 

BLAST 732: Zest = 2.97 (or Zest = 2.63) 

This source does not have a counterpart at 350 /xm or 500 /im, but 
the 250 //m flux-density is too faint for this to offer a useful redshift 
constraint. However, the lack of any radio detection suggests the 
source, if real, lies at high redshift, and indeed the 250 /im position 
lies only 10 arcsec from the fourth brightest 870 /im LAB OCA 
source in GOODS-South (Weiss et al. 2009). We therefore fitted 
galaxy models to the optical-infrared photometry of all objects in 
the ISAAC K-bmd within 15 arcsec of the 250 /xm position. This 
search produced six potential counterparts, three of which have 
estimated redshifts z > 1.5. Of these 3, only one has a 24 /im 
counterpart, so we adopt this as the primary identification. This 
object, which we designate 732-1, is a very red galaxy which 
was too faint to be included in the GOODS -MUSIC catalogue, 
but which fortuitously lies with the Hubble Ultra Deep Field 
recently imaged with WFC3 on HST (McLure et al. 2010). The 
resulting high-accuracy photometry for this source, tabulated in 
Table Al, yields a very robust photometric redshift Zest = 2.97), 
as ilustrated in Fig. 5. The alternative IDs are 732-2, which is 
GOODS MUSIC 30080 (zest = 2.63), and 732-3= GOODS 
MUSIC 10787 (zest = 2.40). Choosing between these alternatives 



clearly does not significantly affect the final redshift distribu- 
tion, and it is certainly possible that both the 250 and 870 /xm 
emission may arise from a blend of a number of sources at 2; c:^ 3.0. 

BLAST 830: Zspec = 0.605 (or Zspec = 0.735) 
This source does not have a counterpart at 350 /im or 500 /im, but 
the 250 /im flux density is too faint for this to offer a useful redshift 
constraint. There are two alternative, statistically- significant radio 
counterparts, and we select VLA 239 in favour of VLA 1179. 
The chosen identification appears to be an interacting galaxy with 
spectroscopic redshift z = 0.605. The edge-on disc galaxy a few 
arcsec NE lies at ^ = 0.735 and so is not physically associated 
with the adopted identification. However, the 250 /im flux density 
may receive contributions from both these objects, and from VLA 
1179. Whatever the exact division of blended emission, it is clear 
that this 250 /im arises primarily from sources at 2; 0.6 — 0.75. 

BLAST 861: Zest = 1 95 

This source does not have a counterpart at 350 /im or 500 /im, but 
the 250 /im flux-density is too faint for this to offer a useful redshift 
constraint. The radio indentification, at a radius of 11.9 arcsec, 
only just passes the adopted significance threshold P < 0.1. This 
identifies the BLAST source with a faint, red, possibly interacting 
galaxy, for which the available photometry delivers a unique and 
well-constrained estimated redshift at 2; c::^ 2. As with BLAST 66, 
the published C0MB017 photometric redshift Zphot = 0.881 is 
very poorly constrained, and can be safely rejected. 

BLAST 983: Zspec = 0.366 

This 250 /im source is not found in the BLAST catalogues at 
longer wavelengths, suggesting it lies at only moderate redshift. 
The radio identication is formally secure with P < 0.1. It lies 
rather far from the 250 /im position, almost at our limiting search 
radius, but this is not unreasonable given that it is the second least 
significant source in our sample. The resulting galaxy identifica- 
tion at 2; = 0.366 can clearly be seen to be dusty in Fig. 7. Only 
3 arcsec distant is another less massive red galaxy at essentially the 
same redshift (z = 0.368). 

BLAST 1293: z^pec = 1 382 

This source does not have a catalogued counterpart at 350 /im 
or 500 /im, but the 250 /im flux-density is too faint for this 
to offer a useful redshift constraint. There are two alternative, 
formally- significant, radio counterparts, and we select VLA 211 
in preference to VLA 148. While the 250 /im source could be 
a blend of both objects, our statistically chosen identification 
transpires to be a bright source at 24 /xm, while the alternative 
radio counterpart is undetected in the mid-infrared. As can be 
seen from Fig. 7, our selected identification is a red, high-redshift, 
apparently interacting galaxy. The GOODS -MUSIC redshift for 
this object is 2; = 1.99, but our own photometric redshift for this 
source is Zphot — 1-37 (see last panel of Fig. 5). The latest release 
of VLT F0RS2 GOODS-South redshifts yields Zspec = 1.382 
(Vanzella et al. 2008). 
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BLAST 109: Zspec = 0.124 






Figure 7. Optical/infared postage stamp images of the primary galaxy counterparts we have identified for the 20 BLAST 250 /im sources in GOODS-South, 
ranked by redshift (with BLAST source number and redshift given below each left-hand stamp). Each row of stamps commences with a 36 x 36 arcsec HST 
ACS z-band image, centred on the BLAST source position, with the identification search area indicated by the 15 arcsec radius circle. White diamonds mark 
the available VLA candidate counterparts, with the selected counterpart marked by the orange circle. The remaining three stamps in each row are centred on 
the position of the selected counterpart, and are 12 x 12 arcsec in size. From left- to-right these show a B + y + 2; HST ACS colour image, the Xs-band 
VLT ISAAC image and the 3.6 /xm Spitzer IRAC image. Display levels on the colour images have been set to +100(7, —la for a. z = source, and then 
these values are reduced by the factor (1 + z)^ at progressively higher redshifts (to offset surface brightness dimming). In all grey-scale plots black is set to 
+8cr, and white to — 8cr where a is the pixel rms. As explained in Section 2, two sources lie just outside the available near-infrared imaging, and hence lack a 
Kg -band postage stamp. 

© 0000 RAS, MNRAS 000, 000-000 



18 J.S. Dunlop, et al. 




BLAST 983: ^spec = 0.366 




BLAST 830-1: Zspec = 0.605 




BLAST 257: ^est = 0.689 



Figure 7. continued 
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BLAST 1293: Zspec = 1.382 




BLAST 552: ^est = 1.68 




BLAST 193: Zest = 1 81 




BLAST 158: ^est = 1.85 



Figure 7. continued 
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Figure 7. continued 
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BLAST 593: ^est > 2.5 



Figure 7. continued 
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BLAST 6-2: Zspec = 0.076 
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Figure 8. Optical/infrared postage stamp images of the alternative galaxy identifications for 6 of the BLAST sources, including two alternative candidates for 
BLAST 732. Individual panels are as described in the caption to Fig. 7. 
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BLAST 732-3: ^est = 2.6 



Figure 8. continued 
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BLAST LABOCA ds-L dL- 



Plb Notes 



ID 


ID 


/ CUCSCV./ 








59 


10+34 


8.9 


7.7 


0.004 


Blend 


66 


18 


7.5 


5.5 


0.008 




158 


79 


7.4 


5.0 


0.011 




318 


67 


12.4 


3.7 


0.019 




593 


12 


8.7 


2.4 


0.008 




732 


32 


9.6 


4.4 


0.012 


Opt ID 



Table 4. Associations between BLAST 250)L^m and LABOCA LESS 
870 /im sources within the central 150 arcmin-^ of GOODS-South. (Ib-l ^^ 
the distance, in arcsec, between the BLAST and LABOCA source. For com- 
parison, dL-Rad is the distance, in arcsec, between the LABOCA source 
and the radio identification listed in Table 2. Plb is the probability that 
the association between the BLAST and LABOCA source is the result of 
chance, calculated in an analogous manner to the P values calculated for the 
BLAST-Radio associations in Table 2, but with a search radius of 18 arcsec. 



6 DISCUSSION 

6.1 LABOCA sub-mm detections and 250/870 /im flux ratios 

Within the central area of GOODS-South under study here, Weiss 
et al. (2009) have extracted 10 LABOCA 870 /im sources with a 
raw significance > 4cr. Searching around the BLAST positions out 
to a search radius of 18 arcsec, we find that seven of these 870 /xm 
sources coincide with 250 //m sources to within a positional accu- 
racy of < 13 arcsec (including two associated with BLAST 59). 
These associations are tabulated in Table 4. We adopted a search 
radius of 18 arcsec by adding the positional uncertainty adopted 
for the BLAST sources in Section 3.1.1 {opos — 6 arcsec) in 
quadrature with an assumed positional uncertainty in the LABOCA 
sources of cFpos = 4 arcsec (calculated assuming the 19 arcsec 
LABOCA beam, and a typical deboosted S/N = 3), and then multi- 
plying by 2.5 to ensure 95% completeness. 

The LABOCA and BLAST source surface densities are both 
so much lower than the radio, mid-infrared or optical source surface 
densities, that all of these associations are statistically compelling, 
as demonstrated by the derived values of Pl b given in Table 4. The 
other 3 LABOCA sources in the field are not near to any BLAST 
source, and so the choice of search radius appears to have been 
sensible. 

The associations listed in Table 4 provide an opportunity to 
check that our adopted positional uncertainties are indeed reason- 
able. First, the median value of dL-Rad = 5.0 implies that (assum- 
ing zero uncertainty in the radio positions), the positional uncer- 
tainty in the LABOCA positions is apos ^4.5 arcsec (allowing for 
the radial weighting in the observed distribution), showing that our 
adoption of 4 arcsec was at least approximately correct. The median 
value of the BLAST-LABOCA positional offset is dB-L = 8.9, 
implying aB-L — 7.5 arcsec. Subtracting 4.5 in quadrature from 
7.5, leads to the conclusion that apos = 6 arcsec for the BLAST 
sources, as adopted in Section 3.1.1. Repeating this analysis with 
means instead of medians, yields apos = 6.5 arcsec. 

Working backwards, this empirical check on the uncertainty 
in the positions of the BLAST sources considered in this paper re- 
affirms that they are indeed genuine :^ 4cr sources prior to flux 
deboosting (including the contribution of confusion to the noise). 

The more scientifically interesting aspect of these LABOCA 
detections is that all of them are associated with high-redshift 
BLAST sources. The 870 )Lxm flux densities of the LABOCA de- 
tections are listed in Table 3. From examination of this Table it 
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Figure 9. The LABOCA LESS 870 /im detections and non-detections of 
the BLAST 250 jim sources. The upper panel shows 870 /im flux den- 
sity plotted against redshift (spectroscopic or photometric) as derived from 
the primary proposed galaxy identification for each BLAST source. All the 
870 /im detections are confined to the high-redshift BLAST subsample. The 
non-detections are indicated by arrows plotted at the adopted conservative 
~ 3-<j flux-density limit of Ssio < 4mJy. The lower panel plots the 
derived 250/870 /xm flux-density ratios (or limits), again versus redshift. 
The anticipated trend of declining 250/870 /im flux-density ratio with red- 
shift is clearly evident. The solid line, derived by redshifting a simple mod- 
ified black-body spectrum with T = 40K, and (5 = 1.5, provides a very 
good description of the data. The dotted curve shows the effect of increas- 
ing temperature to T = 45 K, and the dashed line the effect of reducing it 
to T = 35 K. In calculating the observed flux ratios, neither the BLAST 
nor the LABOCA flux densities have been deboosted. 
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can be seen that the LAB OCA detections are confined to proposed 
identifications with z > 1.5. This point is illustrated in the up- 
per panel of Fig. 9, where we have also adopted a conservative 
~ 3cr flux-density limit of Sgro < 4 mJy for the non-detections. 
Not surprisingly, it can be seen that the redshift distribution of the 
LABOCA-detected BLAST sources is consistent with that which 
is displayed by galaxy samples selected at 850//m with SCUBA. 

The 870//m detections thus provide strong, independent sup- 
port for most of our proposed high-redshift BLAST identifications, 
especially since the identifications and redshifts were established 
without any reference to the LAB OCA data. However, given the 
existence of the LAB OCA data, one can then also ask whether the 
LABOCA wo/i-detections of 4 of our proposed z > 1.5 identifi- 
cations cast doubt on their reality. To check this we have therefore 
plotted, in the lower panel of Fig. 9, the 250/870 //m flux-ratios 
of the BLAST galaxies versus redshift. Here, despite the fact that 
only lower limits are available at low redshift, the anticipated trend 
of declining 250/870 /xm flux-density ratio with redshift is evident. 
It is also clear that, with one exception (BLAST 654), the non de- 
tection of the high-redshift candidates is not sufficiently deep to 
cast serious doubt on the proposed redshifts. 

The solid curve shown in the lower panel of Fig. 9 is the pre- 
dicted decline in colour with increasing redshift for a simple mod- 
ified black-body spectrum with an assumed rest-frame temperature 
T = 40K and dust emissivity emissivity index /3 = 1.5. This of- 
fers a remarkably good description of the high-redshift data and is 
at least consistent with the lower redshift limits. In other words, 
in retrospect we should not be surprised that the 250 /im-selected 
galaxies ai z < 1 have not been detected in the LABOCA LESS 
map, nor that the deep 250 /xm imaging analysed here has reached 
sufficient depth to allow us to detect the sub-mm galaxy popula- 
tion at z ^ 2. The dashed and dotted curves simply illustrate the 
effect of increasing or decreasing the assumed dust temperature by 
5K, keeping /3 = 1.5. Given the mismatch between the BLAST 
and LABOCA beams, and the fact that neither set of flux densi- 
ties has been deboosted in constructing Fig. 9, we do not want to 
over-emphasize the limited extent to which one can determine the 
temperature of the dust emission from the 250/870 /im flux-ratios 
of the BLAST galaxies. However, at the very least. Fig. 9 provides 
a basic sanity check on the plausibility of our inferred galaxy iden- 
tifications and redshifts. 



6.2 The redshift distribution for the 250 fim sample 

The inferred redshift distribution of our 20-source BLAST 250 fim 
sample is presented in Fig. 10. We show two alternative versions, 
in order to illustrate the extent to which its basic form is robust 
against changes in ambiguous identifications. BLAST 654 and 552 
are excluded from both plots, the former because its 250/870 //m 
limit appears anomalously high for its proposed redshift, the latter 
because its non detection with LABOCA, combined with the low 
significance of its radio identification and a non detection at 24 /am 
casts serious doubt on the robustness of this identification. In the 
upper panel we plot the histogram derived from the primary galaxy 
identifications/redshifts listed in Table 3 (we include both 59-1 and 
59-2). In this case 10/19 sources lie dX z > 1, and the distribution 
appears almost bi-modal. In the lower panel we have replaced the 
primary identifications with the secondary/alternative galaxy ID for 
all potentially ambiguous cases (i.e. BLAST sources 6, 257, 503, 
732, 830) and have also simply ejected BLAST 193 from the sam- 
ple (because it possesses neither a VLA radio identification, nor a 
LABOCA sub-millimetre detection). This more sceptical approach 
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Figure 10. An exploration of the robustness of the basic form of the red- 
shift distribution for the 250 jivn sources in GOODS-South. Sources 654 
and 552 have been rejected from both distributions for the reasons detailed 
in the text. The upper panel shows the redshift distribution derived from the 
primary identifications (plus 59-2) listed in Table 3. In the lower panel we 
have replaced the primary identifications by the secondary identifications 
in the 5 appropriate cases, and in addition have removed the other source 
which lacks either a radio or a LABOCA counterpart (BLAST 193). Despite 
these 6 changes, the two distributions are statistically indistinguishable. In 
the upper panel there are 10/19 galaxies at z > 1. In the lower plot there 
are 8/18 sources at 2; > 1. A median redshift 2; ~ 1 is clearly a fair descrip- 
tion of the data presented here, given the uncertainties in identifications and 
redshifts. 
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Figure 11. The location of the adopted primary (red circles) and secondary (green stars) galaxy counterparts for the BLAST 250 /xm sources on the Spitzer 
IRAC colour-colour plane. The small grey points indicate the positions of field galaxies within GOODS-South, as derived from all galaxies in the GOODS- 
MUSIC catalogue down to limiting IRAC AB magnitudes of ms.e = 24.0, 777,4.5 = 23.4, ms.s = 22.0, mg.o = 22.0. Each BLAST identification is 
labelled with its redshift. The three tracks are derived by simulating the effect of redshifting, from z = to z = 3.5, the spectral energy distributions of 
two well-known local star-forming galaxies (the Virgo Sc galaxy VCC1972, and M82) and the more highly obscured ULIRG Apr220. The blue solid track 
is for VCC1972, the green dashed track is for M82, and the red dot-dash track is for Arp220, with the SEDs for all three galaxies taken from the multi- 
wavelength fits derived by Devriendt et al. (1999). The locations of the BLAST primary identifications are in general very well described by the VCC1972 
or M82 tracks, which intially move rapidly in observed 5.8 — 8.0 /im colour as the S/xm PAH feature moves out of the longest IRAC filter, and later move 
rapidly up in 3.6 — 4.5 /im colour as the 1.6/im peak in the starlight moves into the 4.5 /im filter. Very few of the BLAST galaxies lie on the Arp220/ULIRG 
track, with the vast majority being far better described by the less reddened star-forming galaxies. At low redshift this tallies well with the appearance of the 
BLAST identifications shown in Fig. 7, which shows that they are generally isolated dusty star-forming discs. Moreover, certainly at 2; < 1, 250 /im selection, 
especially over the relatively small field of GOODS-South will be biased towards cooler SEDs which peak longwards of A = 100 /im. This is true for 
star-forming disc galaxies, but not true for the hotter ULIRGs selected at 60 /im with IRAS. Perhaps more surprising (given their 'ULIRG-like' star-formation 
rates) is the fact that the z ~ 2 BLAST/LABOCA galaxies also appear relatively unobscured at IRAC wavelengths. 



still leaves 8/18 sources at 2; > 1, of which 6 have been detected by 
LABOCA. 

For our sample, these two alternative redshift distributions 
are statistically indistinguishable, and we conclude that, at least 
in GOODS-South, ^ 50% of the 250 /xm sources with S250 > 
35 mJy lie at 2; > 1. Any evidence for bimodality is statistically 
insignificant, but the median redshift, and the fact that 1 /3 of 
the sample lies at z > 1.5, appear to be robust results. The precise 
form of the distribution is of course aflicted here both by small- 
number statistics, and by the potential effects of cosmic variance. 
It will therefore be interesting to see the extent to which the red- 
shift distribution derived here applies to the much larger samples 
of 250 /im sources which should be uncovered with Herschel. 

It is interesting to assess how our results on this small, deep. 



near-complete BLAST subsample compare with those reported for 
the full 5-cr BLAST sample of the wider 9 degree^ of the Extended 
Chandra Deep Field South by Dye et al. (2009). This comparison 
must be done with care, and the limitations understood, because 
Dye et al. were studying a larger but brighter sample, with on av- 
erage much poorer supporting data. Dye et al. attempted both radio 
and 24 /im identifications for 350 BLAST sources, secured identi- 
fications for 198 of them, and in the end secured redshifts for 74 
of these (i.e. for ~l/5th of the sample). It is not clear what conclu- 
sions can be drawn from the resulting redshift distribution, given 
the low-level of completeness and the mix of flux-density limits. 

More instructive for the present comparison is the redshift dis- 
tribution derived by Dye et al. for the > 5cr BLAST sources in the 
central region of the ECDFS, with the best supporting data from 
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the FIDEL survey. Here Dye et al. identified 50% of the BLAST 
sources (c^ 75% of the 250 )Ltm-selected sources) at either radio 
or mid-infrared wavelengths. The median redshift reported in this 
sub-area by Dye et al. is Zmed — 0.95. Given that the redshift in- 
formation includes some COMBO 17 redshifts which (as discussed 
here for, e.g. BLAST 66) are undoubtedly under-estimates, given 
that the Dye et al. redshift distribution is more incomplete than that 
presented here, and given that the effective 250 //m flux density 
limit of the Dye et al. deep sample is still > 50% higher than in the 
present study, we conclude there is no real evidence for any incon- 
sistency between our redshift distribution and that derived by Dye 
et al. (2009). The only visible difference is that the high-redshift 
tail in the present deep study extends to somewhat higher redshifts 
than that reported by Dye et al., but this is unsurprising given the 
above-mentioned selection effects. 

6.3 The mid-infrared colours of the 250 fim galaxies. 

Finally, in Fig. 11 we briefly investigate the location of the galaxy 
counterparts for the BLAST 250 /xm sources on the Spitzer IRAC 
colour-colour plane. This is of interest for a number of reasons. 
First, it is worthwhile to check how the mid-infrared colours of 
the BLAST-selected galaxies compare with those of the general 
galaxy population in GOODS-South. Second, we can compare their 
colours with expectations based on the redshifted SEDs of known 
local far-infrared bright galaxies. Third, it is of potential interest to 
check whether the counterparts of 250 /im-selected galaxies can, 
in future, be reliably isolated on the basis of IRAC colour (a point 
of importance given the potentially extensive sky coverage which 
could be provided by the Spitzer warm mission). 

We choose to plot 3.6 iim. — 4.5 /xm colour versus 5.8 //m — 
8.0 //m colour, as advocated by Stern et al. (2005), rather than 
3.6 jim. — 5.8 jim. colour versus 4.5 /xm — 8.0 jim. colour as ad- 
vocated by some other authors (e.g. Lacy et al. 2004; Pope et al. 
2006). We do this to avoid unhelpful correlations, and because only 
3.6 /im — 4.5 /im colours will be available from the Spitzer warm 
mission. 

In Fig. 1 1 we plot the location of the adopted primary and sec- 
ondary galaxy counterparts for the BLAST 250 jim. sources, super- 
imposed on the positions of field galaxies within GOODS-South, 
as derived from all galaxies in the GOODS -MUSIC catalogue 
down to limiting IRAC AB magnitudes of ms.e = 24.0, m4.5 = 
23.4, 7715.8 = 22.0, ms.o = 22.0. In addition, we show tracks in 
the colour-colour diagram produced by redshifting (from z = to 
z = 3.5) the spectral energy distributions of two well-known local 
star-forming galaxies (the Virgo Sc galaxy VCC1972, and the star- 
burst galaxy M82), along with the analogous predicted track for the 
highly-obscured ULIRG Apr 220. 

The main points we can learn from this diagram can be sum- 
marised as follows. First, the primary BLAST galaxy identifica- 
tions are very well described by the two star-forming galaxy SEDs, 
with the red data points essentially mapping out the redshift evo- 
lution of the tracks. The rapid horizontal movement at low redshift 
(z = — 0.5) is caused by the poly cyclic aromatic hydrocarbon 
(R\H) features moving through the longer two IRAC bands. The 
subsequent rapid vertical movement from 3.6 fim — 4.5 /xm < to 
3.6 /im — 4.5 //m > is the result of the stellar photospheric peak 
at Xrest — 1.6//m moving through the shorter two IRAC bands at 
z 1.5. 

Second, it can be seen that very few of our primary iden- 
tifications appear to be analogues of Arp220, although at the 
very highest redshifts all three templates become indistinguish- 



able. Our low-redshift BLAST galaxies stay relatively blue in 
3.6 jjLin — 4.5 jjLin colour, indicating their star-light is not swamped 
by the mid-infrared power-law produced by very hot dust (AGN or 
starburst heated), as is seen in Arp 220 and other local ULIRGs. 
Even 2Li z 2 most of the BLAST/LABOCA galaxies have 
colours expected from starlight (albeit in some cases reddened to 
3.6 fiui — 4.5 fim > 0.25). Our highest-redshift candidates do lie 
in the region of the diagram expected for ^ :^ 3, but at this point 
essentially all tracks converge on the same location in the diagram. 

It is interesting to speculate why the Arp 220 track is such 
a poor description of the locus occupied by the BLAST sources. 
At low redshift it is perhaps unsurprising, as starburst galaxies 
like M82 and Sc star-forming discs both have (relatively) cool 
SEDs which peak longward of Xrest — 100 /xm, whereas the 
ultra-luminous ULIRGS first uncovered by IRAS at 60 /im peak at 
shorter wavelengths. A 250 jum survey is thus biased towards the 
selection of these cooler more 'normal' starbursts at low redshift, 
especially since highly-obscured ULIRGs are rather rare objects. 
The appearance of the low-redshift BLAST galaxies in Fig. 7 (i.e. 
star-forming discs) is thus consistent with their IRAC colours. 

Perhaps more surprising (given their 'ULIRG-like' star- 
formation rates) is the fact that the 2; :^ 2 BLAST/LABOCA galax- 
ies also appear relatively unobscured at IRAC wavelengths. How- 
ever, this is arguably consistent with their 250/870 //m flux-ratios 
as discussed in Section 6.1, which appear to be consistent with a 
relatively cool dust temperature of 40 K (see Fig. 9). 

Finally we revisit the claim made by Yun et al. (2008) that a 
simple and robust sub-mm galaxy candidate selection criterion is 
provided by 3.6 /im — 4.5 /im > —0.2. It is obvious from Fig. 11 
that this criterion is in fact satisfied by essentially any galaxy at 
z > 1.5, and thus cannot be used to isolate any special sub-class of 
source at these redshifts. 

These issues, along with the other physical characteristics of 
the BLAST galaxies (size, morphology, stellar mass, star-formation 
history) will be explored in more detail by Targett et al. (in prepa- 
ration). 



7 SUMMARY 

We have identified and investigated the nature of the 20 bright- 
est 250 /im sources detected by BLAST within the central 150 
arcmin^ of the GOODS-South field. Aided by the available deep 
VLA imaging, reaching Si a — 40/iJy (4-cr), we have identified 
secure radio counterparts for 17/20 of the 250 /xm sources. The re- 
sulting enhanced positional accuracy of c± 1 arcsec has then al- 
lowed us to exploit the deep optical (HST), near-infrared (VLT) 
and mid-infrared (Spitzer) imaging of GOODS-South to establish 
secure galaxy counterparts for the 17 radio-identified sources, and 
plausible galaxy candidates for the 3 radio-unidentified sources. 

Confusion is a serious issue for surveys with such large beams, 
and in many cases there may be no such thing as a unique, single 
counterpart. We can nevertheless expect our chosen counterparts to 
be significant, and often dominant contributors to the BLAST flux 
densities, and we give several arguments to support this claim. It 
is therefore still reasonable to pursue the properties of this 'bright- 
est counterpart' sample. For all of our 20 chosen identifications we 
have been able to determine spectroscopic (8) or robust photomet- 
ric (12) redshifts. The result is the first near-complete redshift dis- 
tribution for a deep 250 /xm-selected galaxy sample. This reveals 
that 250 fxm surveys reaching detection limits S250i^ui — 40 mJy 
contain not only low-redshift spirals/LIRGs, but also the extreme 
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z :^ 2 dust-enshrouded starburst galaxies previously discovered 
at longer sub-millimetre wavelengths. Inspection of the LAB OCA 
870 /xm imaging of the GOODS-South field yields detections of 7 
of the proposed z > 1 BLAST sources, and reveals 250/870 iim 
flux ratios consistent with a standard 40K modified black-body fit 
with a dust emissivity index /3 = 1.5. Thus, at least in GOODS- 
South, we infer that ^ 50% of the 250 /im sources with 5'25o > 
40mJy lieatz > 1. 

Based on their IRAC colours, we find that virtually all of 
the BLAST galaxy identifications appear better described as ana- 
logues of the M82 starburst galaxy, or Sc star-forming discs rather 
than highly obscured ULIRGs. This is perhaps as expected at low- 
redshift, where the 250 /xm BLAST selection function is biased to- 
wards SEDs which peak longward of Xrest — 100 /xm. However, 
it also appears largely true at z 2. 



ACKNOWLEDGEMENTS 

JSD acknowledges the support of the Royal Society through a 
Wolfson Research Merit Award, and the support of the European 
Research Concil through the award of an Advanced Grant. We ac- 
knowledge the support of NASA through grant numbers NAG5- 
12785, NAG5-13301, and NNG0-6GI11G, the NSF Office of Po- 
lar Programs, the Canadian Space Agency, the Natural Sciences 
and Engineering Research Council (NSERC) of Canada, and the 
UK Science and Technology Facilities Council (STFC). This work 
is based in part on observations made with the Spitzer Space 
Telescope, which is operated by the Jet Propulsion Laboratory, 
California Institute of Technology under a contract with NASA. 
APEX is operated by the Max-Planck-Institut fur Radioastronomie, 
the European Southern Observatory, and the Onsala Space Ob- 
servatory. This work is based in part on observations made with 
the NASA/ESA Hubble Space Telescope, obtained from the Data 
Archive at the Space Telescope Science Institutewhich is operated 
by the Association of Universities for Research in Astronomy, Inc., 
under NASA contract NAS 5-26555. MC acknowledges the award 
of a STFC Advanced Fellowship. IRS acknowledges support from 
STFC. JLW acknowledges the support of an STFC Studentship. 



REFERENCES 

Afonso J., et al., 2006, ApJ, 131, 1216 

Austermann J.E., et al., 2010, MNRAS, 401. 160 

Beckwith S.V.W., et al., 2006, AJ, 132, 1729 

Bertoldi P., et al., 2007, ApJS, 172, 132 

Bouwens R.J., et al, 2008, ApJ, 686, 230 

Bouwens R.J., et al., 2010, ApJ, 709, L133 

Bolzonella M., Miralles J.-M., Pello R., 2000, A&A, 363, 476 

Brammer G.B., van Dokkum RG., Coppi R, 2008, ApJ, 686, 1503 

Bruzual G., 2007, ASPC, 374, 303 

Calzetti D., et al, 2000, ApJ, 533, 682 

Caputi K.L, McLure R.J., Dunlop J.S., Cirasuolo M., Schael A.M., 2006, 

MNRAS, 366, 609 
Chapin E.L., et al., 2010, MNRAS, submitted <arXiv:1003.26?7l 
Coppin K., et al., 2006, MNRAS, 372, 1621 
Devlin M.J., et al., 2009, Nature, 458, 737 

Devriendt J.E.G., Guiderdoni, B., Sadat, R., 1999, A&A, 350, 381 
Dickinson M., et al., 2004, ApJ, 600, 99 
Doherty M., et al, 2005, MNRAS, 361, 525 

Downes A.J.B., Peacock J.A., Savage A., Carrie D.R., 1986, MNRAS, 218, 
31 

Dunlop J.S., Cirasuolo M., McLure R.J., 2007, MNRAS, 376, 1054 



Dunlop J.S., et al., 1989, MNRAS, 238, 1171 
Dye S., et al., 2009, ApJ, 703, 285 
Bales S., et al, 2009, ApJ, 707, 1779 
Fazio GG, et al., 2004, ApJS, 154, 10 

Fixsen D.J., Dwek E E., Mather J.C., Bennett C.L., Shafer R.A., 1998, ApJ, 
508, 123 

Giavalisco, M., et al., 2004, ApJ, 600, L93 

Grazian A., et al., 2006, A&A, 449, 951 

Griffin M., et al., 2007, Advances in Space Research, 40, 612 

Gusten R., et al, 2006, A&A, 454, 13 

Holland W., et al., 2006, SPIE, 6275, 45 

Hughes D.H., et al., 1998, Nature, 394, 241 

Ibar E., Ivison R.J., Biggs A.D., Lai D.V., Best RN., Green D.A., 2009, 

MNRAS, 397, 281 
Ivison R.J., et al, 2002, MNRAS, 337, 1 
Ivison R.J., et al., 2007, MNRAS, 380, 199 
Ivison R.J., et al., 2010, MNRAS, 402, 245 
Kellermann K., et al., 2008, ApJS, 179, 71 
Lacy M., et al., 2004, ApJS, 154, 166 
Le Fevre O., et al., 2004, A&A, 428, 1043 
Le Fevre O., et al, 2005, A&A, 439, 845 
MadauR, 1995, ApJ, 441, 18 
Magnelli B., et al., 2009, A&A, 496, 57 
Marsden G, et al., 2009, ApJ, 707, 1729 

McLure R.J., Dunlop J.S., Cirasuolo M., Koekemoer A.M., Sabbi E., Stark 

D.R, Targett T., Ellis R.S., 2010, MNRAS, 403, 960 
Mignoli M., et al, 2005, A&A, 437, 883 
Miller N.A., et al, 2008, ApJS, 179, 71 
Norris R.R, et al., 2006, AJ, 132, 2409 
Oesch RA., et al., 2010, ApJ, 709, L21 
Pascale E., et al., 2008, ApJ, 681, 400 
Pascale E., et al., 2009, ApJ, 707, 1740 
Patanchon G, et al., 2009, ApJ, 707, 1750 
Peacock J.A., et al., MNRAS, 2000, 318, 535 
Pope A., et al., 2006, MNRAS, 370, 1185 

Puget J.-L., Abergel A., Bernard J.-P, Boulanger F, Burton W.B., Desert 

F-X., Hartmann D., 1999, A&A, 308, L5 
Ravikumar CD., et al., 2007, A&A, 465, 1099 
Retzlaff J., et al, 2010, A&A, 511, 50 

Roche N.D., Dunlop, J.S., Caputi K.L, McLure R.J., Willott C.J., Crampton 

D., 2006, MNRAS, 370, 74 
Rowan-Robinson M., et al., 2008, MNRAS, 386, 697 
Scott K.S., et al., 2010, MNRAS, 405, 2260 
Scott S., et al., 2002, MNRAS, 331, 817 
Serjeant S., et al. 2003, MNRAS, 344, 887 
Stanway E.R., et al., 2004, ApJ, 607, 704 
Stern D., et al., 2005, ApJ, 631, 163 
Strolger L.-G, et al., 2004, ApJ, 613, 200 
Szokoly GP, et al., 2004, ApJS, 155, 271 
Taylor E.N., et al., 2009, ApJS, 183, 295 
Thompson R.L, et al., 2005, AJ, 130, 1 
Truch M.D.P, et al., 2009, ApJ, 707, 1723 
van der Wei A., et al., 2004, ApJ, 601, 5 
Vanzella E., et al., 2005, A&A, 434, 53 
Vanzella E., et al., 2006, A&A, 454, 423 
Vanzella E., et al., 2008, A&A, 478, 83 
Vanzella E., et al., 2009, ApJ, 695, 1163 
Viero M., et al., 2009, ApJ, 707, 1766 
Weiss A., et al., 2009, ApJ, 707, 1201 
Wolf C, et al., 2004, A&A, 421, 913. 

Wolf C, Hildebrandt H., Taylor E.N., Meisenheimer K., 2008, A&A, 492, 
933 

Wuyts S., et al., 2008, ApJ, 682, 985 
Yun M.S., et al., 2008, MNRAS, 389, 333 



© 0000 RAS, MNRAS 000, 000-000 



BLAST 250 jim galaxies in GOODS-South 

APPENDIX A: OPTICAL-INFRARED PHOTOMETRY 

Table Al provides the optical-infrared photometry which was used 
to derive the BLAST galaxy photometric redshifts. 
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BLAST 
ID 


RA 

(J2000) 


Dec 
(J2000) 


435 


606 


775 


^850 


J 


H 


K 


3.6 


4.5 


6 


8.0 


'S'24 




4 


53.146168 


-27.925831 


17.59±0.01 


16.57±0.01 


15.97±0.01 


15.62±0.01 


15.11±0.01 




14.83±0.01 


15.41±0.01 


15.88±0.01 


15.38±0.01 


14.11±0.01 


7090±709 


6-1 


53.124493 


-27.740086 


17.57±0.01 


16.86±0.01 


16.46±0.01 


16.30±0.01 


16.01±0.01 


15.87±0.01 


15.92±0.01 


16.25±0.01 


16.65±0.01 


16.22±0.01 


14.16±0.01 


7590±759 


109 


53.074467 


-27.849800 


20.36=b0.01 


19.10±0.01 


18.41±0.01 


18.07±0.01 


17.48±0.01 


17.11±0.01 


17.00±0.01 


17.43±0.01 


17.75±0.01 


17.83±0.01 


15.95±0.01 


1270±127 


637-1 


53.183559 


-27.862030 


21.11±0.01 


19.61±0.01 


18.82±0.01 


18.46±0.01 


17.59±0.01 




16.78±0.01 


17.65±0.01 


17.70±0.01 


18.21±0.01 


17.09±0.01 


594± 36 


983 


53.136730 


-27.768833 


22.49±0.02 


21.10±0.01 


20.39±0.01 


20.03±0.01 


19.41±0.01 


18.86±0.01 


18.30±0.01 


18.87±0.01 


18.77±0.01 


19.13±0.01 


18.02±0.01 


611± 


10 


104 


53.150726 


-27.825510 


22.77±0.01 


21.49±0.01 


20.62±0.01 


20.24±0.01 


19.67±0.01 


19.06±0.01 


18.92±0.01 


19.09±0.01 


19.25±0.01 


19.29±0.01 


19.10±0.01 


744 ± 


7 


830-1 


53.055141 


-27.711374 


21.83±0.01 


20.92±0.01 


20.06±0.01 


19.84±0.01 


19.25±0.01 




18.63±0.01 


18.80±0.01 


18.99±0.01 


18.87±0.01 


19.13±0.02 


774 ± 


10 


257 


53.075142 


-27.831411 


23.23zb0.02 


21.68±0.01 


20.49±0.01 


20.06±0.01 


19.43±0.01 


18.98±0.01 


18.53±0.01 


18.71±0.01 


19.10±0.01 


19.23±0.01 


19.84±0.01 


394± 


5 


1293 


53.193047 


-27.890892 


26.44it0.20 


26.11±0.11 


25.03±0.08 


24.27±0.05 


22.98±0.09 




21.45±0.06 


20.07±0.05 


20.10±0.05 


20.00±0.07 


20.29±0.10 


187± 


4 


552 


53.052277 


-27.718325 


28.34±0.54 


27.39±0.21 


26.23±0.14 


25.74±0.11 


24.15±0.09 




22.84±0.06 


22.09±0.08 


21.66±0.11 


22.07±0.92 


22.16±1.09 


290± 


4 


193 


53.053574 


-27.778013 


28.13zb0.42 


28.70±0.65 


27.18±0.31 


26.27±0.17 


23.98±0.13 


23.29±0.07 


22.58±0.06 


21.53±0.04 


21.32±0.08 


2 1.25 ±0.42 


21.95±0.98 






158 


53.089821 


-27.940001 


25.73±0.09 


24.85 ±0.05 


24.05 ±0.04 


23.59±0.04 








19.83±0.04 


19.55±0.04 


19.64±0.07 


20.17±0.13 


586± 


6 


66 


53.020180 


-27.779888 


27.73zb0.18 


27.13±0.11 


26.56±0.15 


25.43±0.11 








20.84±0.11 


20.35±0.11 


20.21±0.11 


20.83±0.12 


620± 


7 


861 


53.198280 


-27.747866 


27.17±0.15 


26.22±0.11 


25.52±0.12 


24.92±0.07 


23.44±0.09 


22.72±0.11 


22.30±0.08 


20.99±0.08 


20.83±0.08 


20.55±0.15 


20.93±0.21 


299 ± 


4 


503-1 


53.157238 


-27.833452 


25.90±0.05 


25.26±0.02 


24.55±0.02 


24.18±0.02 


23. 15 ±0.04 


22.46±0.03 


22.13±0.03 


21.07±0.02 


20.73 ±0.02 


20.80±0.05 


21.21±0.08 


251± 


5 


318 


53.179874 


-27.920731 


25.06±0.06 


24.54±0.04 


24.05 ±0.05 


23.80±0.04 


22.90±0.08 




21. 63 ±0.05 


20.56±0.04 


20.09 ±0.05 


19.90±0.05 


20.28 ±0.05 


554± 


5 


59-1 


53.079300 


-27.870552 


25.90±0.08 


25.55±0.06 


24.94±0.06 


24.75±0.07 


24.48±0.15 




23.78±0.14 


23.29±0.11 


22.45 ±0.10 


22.48±0.30 


22.30±0.29 


148± 


5 


654 


53.136547 


-27.885691 


25.41 ±0.07 


24.51 ±0.03 


24.29±0.05 


24.21 ±0.05 


23.81±0.08 




23.10±0.09 


22.73±0.06 


22.53±0.08 


22.50±0.23 


22.82±0.32 






732-1 


53.181458 


-27.777472 


29.94±0.50 


28.23±0.12 


27.98±0.10 


27.60±0.15 


26.54±0.07 


25.36±0.06 


24.04±0.18 


22.25±0.12 


21.78±0.14 


21.44±0.21 


21.16±0.20 


85± 


4 


593 


53.199965 


-27.904560 


> 27.3 


> 27.3 


> 26.5 


> 26.2 


> 25.5 


> 24.7 


24.10±0.20 


21.93±0.06 


21.71±0.13 


21.53±0.22 


21.19±0.22 


49 ± 


3 


6-2 


53.124950 


-27.734684 


18.16±0.01 


17.39±0.01 


16.95±0.01 


16.76±0.01 


16.45±0.01 


16.28±0.01 


16.32±0.01 


16.73±0.01 


17.13±0.01 


16.86±0.01 


14.95±0.01 


4660±466 


503-2 


53.161751 


-27.832291 


21.99±0.01 


20.69±0.01 


19.96±0.01 


19.64±0.01 


18.90±0.01 


18.35±0.01 


17.93±0.01 


18.87±0.01 


18.77±0.01 


19.23±0.01 


17.55±0.01 


543 ± 


5 


637-2 


53.184448 


-27.861420 


21.64±0.01 


20.19±0.01 


19.47±0.01 


19.15±0.01 


18.48±0.01 




17.92±0.01 


18.44±0.01 


18.17±0.01 


17.96±0.01 


17.06±0.01 






830-2 


53.057747 


-27.713591 


23.41 ±0.02 


22.36±0.01 


21.27±0.01 


20.83±0.01 


20.04±0.01 




19.22±0.01 


19.56±0.01 


19.82±0.01 


19.88±0.02 


19.92±0.02 


333± 


7 


59-2 


53.071545 


-27.872448 


24.97±0.06 


23.87±0.02 


22.65±0.01 


21.80±0.01 


21.03±0.01 




19.93±0.01 


19.68±0.01 


19.66±0.01 


19.97±0.02 


19.79±0.02 


288 ± 


6 


732-2 


53.180542 


-27.779686 


28.22±1.09 


27.07±0.17 


27.81±0.72 


27.32±0.52 


26.08±0.54 


24.88±0.30 


23.69±0.09 


22.70±0.08 


22.39±0.16 


21.90±0.51 


21.80±0.65 






732-3 


53.182018 


-27.779537 


25.34±0.07 


24.72±0.03 


24.51 ±0.06 


24.39±0.05 


23.95±0.08 


23.36±0.08 


23.30±0.07 


23.29±0.07 


23.16±0.09 


23.01 ±0.26 


24.19±0.73 







Table Al. Optical (H ST ACS, near-infrared (VLT ISAAC) and mid-infrared (Spitzer IRAC + MIPS) photometry for the BLAST galaxy counterparts. The table lists the primary identifications, ranked by redshift (as 
shown in Figure 7), and then the secondary identifications (as shown in Figure 8). All measurements are given in AB magnitudes, except for the lAjim. photometry, which is given in micro- Jy (the lAfim information 
was not used in the derivation of photometric redshifts, and is simply given here for completeness). Non-detections are given as 2<j limiting magnitudes in 2-arcsec diameter apertures. Where no entry is given, this 
means that the available deep imaging coverage at that particular wavelength did not extend to include the object in question. Note that the optical and J and H photometry for 732-1 comes from the ultra-deep 
imaging of the HUDF provided by H ST ACS and WFC3 imaging (McLure et al. 2010). 



